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2 ' Abstract. We consider the problem of parameter estimation for an ergodic diffusion with Fisher-Snedecor 

invariant distribution, to be called Fisher-Snedecor diffusion. We compute the spectral representation of 
its transition density, which involves a finite number of discrete eigenfunctions (Fisher-Snedecor polyno- 
^ . mials) as well as a continuous part. We propose moments based estimators (related to the Fisher-Snedecor 

Q^ . polynomials) and prove their consistency and asymptotic normality. Furthermore, we propose a statisti- 

O^ . cal test for the distributional assumptions on the marginal distribution of the Fisher-Snedecor diffusion, 

~~^ , based on the moment condition derived from the corresponding Stein's equation. 
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1 Introduction 

Historical roots. The study of diffusion processes with invariant distributions from the Pearson fam- 
ily started in the 1930's, when Kolmogorov |44l [67] studied the Fokker-Planck or forward Kolmogorov 
equation 

-jT = -Tj-{{aix + ao)p) + J— {{b2x'^ + bix + bo)p) , p=p{x,t), x £ R, i > 0, 

with a linear drift and a quadratic squared diffusion, and observed that the invariant density p(-) satisfies 
the differential equation 

p'{x) {ai - 2b2)x + {aa ~ hi) cix + Cq q{x) 



p{x) b2x'^ + bix + bo b2x'^ + bix -\- bo s{x) 



xeR, (1) 



^Corresponding author. 



introduced by K. Pearson [S5] in 1914 (and of illustrious history - see Diaconis and Zabel [S^), in order 

to unify some of the most important statistical distributions. 

It seems appropriate to call this important class of processes Kolmogorov-Pearson (KP) diffusions, or 

Gauss-hypergeometric diffusions, due to the appearance of the Gauss 2^1 function (and its limiting 

confluent forms) in various explicit formulas. 

For a long period of time after that, KP diffusions were neglected, with some notable exceptions like 

Wong (1964) [71], who reemphasized the importance of this class of models as a most natural extension 

of the "first order statistical description characterized by p(a;)" to a time dependent model, and computed 

spectral representations of the transition density in some cases (but missed the Fisher- Snedecor case). 

Mathematical finance motivations. Recently, the interest in these processes was reawakened in the 
context of financial modeling. The most famous case is the Merton-Black-Scholes SDK 

dXt = rXt dt + aXt dWt ^ Xt ^ Xae^"-"' ^^'>'+''^\ (2) 

This has had a huge impact in mathematical finance due to its tractability, which lead to a large variety 
of explicit formulas for the transition and first passage probabilities necessary for the pricing of options. 
Other tractable diffusions are the Ornstein-Uhlenbeck (Vasicek), and Cox-Ingersoll-Ross (CIR) models. 
In the search for more flexible models, financial mathematicians have rediscovered recently the con- 
venience of the Kolmogorov-Pearson diffusions (and others, like diffusions with constant elasticity of 
variance). Some notable contributions are due to Albanese, Campolieti, Carr and Lipton [8], Linetsky 
[5^ [5?1 133] , Kuznetsov [5] 27] , Mendoza, Carr and Linetsky [31] and Shaw [66 . In parallel, the statis- 
tical analysis of these processes was developed by Kiichler and S0rensen and Forman and S0rensen [32] . 
Recently, the interest in developing tractable processes resulted in the introduction by Cuchiero, Keller- 
Ressel and Teichmann [26) of the unifying family of polynomial Markovian processes, whose generator 
maps polynomials into polynomials of (at most) the same degree. 

The Student parametrization of Kolmogorov-Pearson (KP) diffusions. Following the modern 
formulations of [66] and [32], consider the SDE 

dXt ^ 9{fi - Xt) dt + ai dVF^^' + ^2 XtdWJ^^ , t > 0, ^ (3) 

dXt ^e{fi-Xt) dt + 



\ 



>-'+^S^<'-''''(S))«' 



where fc = cr|, Wj and W^ , t > 0, are standard Brownian motions with correlation p, and {Wt, t > 0} 

is a standard Brownian motion resulting from combining the two. 

Notes: 1) The second formulation, to be called the skew Student parametrization, may be used 

for all KP diffusions (by restricting if necessary to the range of values of x where the square root makes 

sense). In the first formulation, one implicitly assumes \p\ < 1,A; > 0, which characterize the Student 

subclass of KP diffusions -see below. 

2) The linearity of the drift and the quadratic variance ensure the existence of polynomial eigenfunctions. 

The S0rensen/skew Student parametrization. Putting p,' = p^,k = -j^, we arrive at 



dXt = 9i^l - Xt) dt + \j —-J {Xt-^if + s^ dBt, i>o, (4) 



where 6^ — (1 — p^)(— )^. This parametrization makes sense for the whole KP faniily (by aUowing 5^ < 0), 
but it will only produce diffusions living on (—00,00) when (5 e M,A; > 0. 



A classification of Kolmogorov-Pearson diffusions in six basic subfamilies may be achieved by 
using the criteria based on the degree deg(s) of the polynomial s{x) from the dcnonrinator of the Pearson 
equation ([1} and on the sign of the leading coefficient k and the discriminant A(s) in the quadratic case. 
The classification is given in the following table 



Pearson diffusion 


Characteristic property 


Invariant density 


Parameter space 


Ornstein-Uhlenbeck 


deg(s) = 


1 

<7%/2 


(x-M)^ 




M e M, 0- > 


-e 2,^2 


Gamma or CIR 


deg(s) = 1 


n,r' '' "^'(o.^)*^' 


a > 0, /3 > 


Beta or Jacobi 


A: < 0, deg(s) = 2, A(s) > 


1 ^_ 


1(1 -x)" ll[o,lj(x) 


a > 0, /3 > 


B(a,3) 


Student (see Q) 


fc > 0, deg(s) = 2, A(s) < 


„.p(ii^.„t«(£^)) 


a>0, 5>0 




h(^)1 


^+1 


Reciprocal gamma 


A: > 0, deg(s) = 2, A(s) = 




a > 0, /3 > 1 


Fisher-Snedecor 


fc > 0, deg(s) = 2, A(s) > 


B(^,|)^' ( = . + « 2 :<„,^j,., 


a > 2, /3 > 2 



Table 1. Classification of ergodic stationary Pearson diffusions. 

Example: In the important Student case, the scale and invariant speed densities are: 

^ ^ arctg (x) 



s(x) = (5(£2 + l))2^e-''^"^*s(*), m{x) 



,x G 



(x2 + l)5k+l 

where we put x := {x — IJ.')/S, 2a — -g — -jzj- 

When 1/ > 1 <^ a > 0, the KP diffusions will have Student stationary distribution with h' degrees 

of freedom, since the speed density may be normalized, arriving thus to diffusions with the invariant 

density: 



f{x) = c{n,fi',a,6)- 
where (see Avram et al. [12] ) 






exp {^ arctg (^)} 



X e 



(5) 



^V^^ (^2 + ^a) tJo 



2aS 



2a 



(6) 



The Ornstein-Uhlenbeck process, CIR process and Jacobi diffusion are well studied and widely applied. 
However, the first results on statistical analysis of reciprocal gamma. Student and Fisher-Snedecor dif- 
fusion, which all have heavy-tailed invariant distributions, are quite new (see [32l |49l [50] ) . Their study 
involves the analysis of the spectrum of the corresponding infinitesimal generators. Namely, the spectrum 
of the infinitesimal generator of Ornstein-Uhlenbeck process, CIR process and Jacobi diffusion is sim- 
ple and purely discrete with classical orthogonal polynomials as corresponding eigenfunctions: Hermite, 
Laguerre and Jacobi polynomials, respectively. In the case of Pearson diffusions with heavy-tailed invari- 
ant distributions, the spectrum of the infinitesimal generator consist of two disjoint parts: the discrete 
spectrum (consisting of finitely many simple eigenvalues) and the essential spectrum. Furthermore, in all 
these cases corresponding eigenfunctions are less known finite systems of orthogonal polynomials: Bessel 



polynomials for reciprocal gamma diffusion (see |49) ) . Routh-Romanovski polynomials for Student diffu- 
sion (see |50] ) and polynomials related to Fisher-Snedecor invariant distribution which have no common 
name (we will refer to them as to Fisher-Snedecor polynomials). 

In this paper we focus on the statistical analysis of the ergodic stationary Fisher-Snedecor diffusion, by 
which we mean solutions of the non-linear stochastic differential equation 



dXt = -e {Xt -fi)dt + J2ext [ ^,5* ^ +J^\dWt, t> 0, (7) 



J/2-1 a/2 
with speed measure/invariant distribution proportional to: 

— — ITTT I(o,oo> (a:) (8) 

(aa;-|-/i(/3-2))2+f 

By taking /^ as the mean of the Fisher-Snedecor distribution pL = a^, (see Section [2]) which may be 
always achieved by scaling, and by assuming a > and P > 2 (ensuring thus the existence of the first 
moment), we arrive to the classic Fisher-Snedecor distribution FS{a,j!i) with (a,/3) degrees of freedom 
as invariant distribution. 

Statistical analysis of stochastic processes often demands the acquaintance of the transition density. 
Since for the most of diffusion models transition density is unavailable in the closed form, the statistical 
analysis is often very complicated. A possible solution to this problem is estimation of the transition 
density, which is used in a recent paper by Beskos [17) in the likelihood based framework. Ai't-Sahalia and 
Mykland [51 E] , in the approach based on the moment conditions which satisfy certain differentiability 
conditions, used approximation of the asymptotic variance by the Taylor expansion with terms known in 
the closed form. However, since for the Fisher-Snedecor diffusion the structure of the spectrum of the 
infinitesimal generator can be fully determined, the use of the spectral representation of the transition 
density is implied. Therefore we present in Section 4 a closed expression for the spectral representation 
of the transition density of Fisher-Snedecor diffusion given in form of the finite sum related to discrete 
spectrum of the infinitesimal generator (i.e. eigenvalues and eigenfunctions - Fisher-Snedecor polynomials) 
and the integral which is taken over the absolutely continuous spectrum of the infinitesimal generator. The 
complicated form of this expression has significant impact on methods used in the statistical part of this 
paper. However, orthogonality of Fisher-Snedecor polynomials and hypergeometric functions appearing 
in the continuous part of the spectral representation makes this result applicable in the statistical analysis 
(see Remark l4.8p . However, in statistical applications treated here in Sections 5 and 6, the discrete part 
of the spectral representation is used much more frequently than the continuous part. This suggests the 
importance of Fisher-Snedecor polynomials which have been ignored in standard mathematical books, 
for instance in Abramowitz and Steghun [Ij, Chichara [27], Erdely [3T], Nikiforov and Uvarov [61] and 
Szego [71] . Significant progress in their research is recently made by Masjed- Jamei [56l [57] and by Koepf 
and Masjed- Jamei [13]. Note that since essentially only the discrete part of the spectrum is used here, 
readers motivated mainly by statistical applications may skip Sections 4.2. - 4.5. 

In statistical part of the paper we observe the problems of parameter estimation and testing hypothesis 
about invariant distribution of the ergodic Fisher-Snedecor diffusion. Note that statistical inference for 
ergodic diffusions is widely studied (see the recent book by Kutoyants [IS], the paper by Kutoyants 
and Yoshida [46] and references therein). The recent review on diffusion estimation can be found in a 
survey paper by H. S0rensen [69]. Significant progress in statistical analysis of Pearson diffusions was 
recently made by Forman and S0rensen |32j , where the problem of parameter estimation is treated by 



method based on the martingale estimation equations. It is proved that the resuhing estimators are 
consistent and asymptoticahy normal with given explicit expressions for the elements of the limiting 
covariance matrix. However, martingale estimation equation does not always provide estimators in the 
explicit form and calculation of the covariance matrix can be quite complicated. Furthermore, in our 
case the maximum likelihood method and the Bayesian estimation theory are not available since the 
diffusion coefficient depends on unknown parameters. Also, the minimum contrast estimation procedure 
based on Whittle functionals (see Anh et al. [TU] and Leonenko and Sakhno |35]) is not applicable here, 
since the fourth order spectral density could not be explicitly calculated. Therefore, parameters of the 
Fisher- Snedecor diffusion are estimated by the simple method of moments. As shown in Kutoyants and 
Yoshida [33] , asymptotic efficiency of the method of moments can be quite good. Another reason that 
makes method of moments suitable for our problem are explicit forms of resulting estimators. Their 
consistency follows from the ergodic property of the observed diffusion. Furthermore, their asymptotic 
normality is implied by the functional central limit theorem based on the strong mixing property with an 
exponentially decaying rate (see Genon-Catalot et al. i34j) and the functional delta method (see Serfling 
[65|). Explicit form of the limiting covariance matrix of the multivariate estimator of parameters of 
invariant distribution is calculated according to the new method based on the closed form expression for 
the spectral representation of the transition density and the finite system of orthogonal Fisher-Snedecor 
polynomials. This method basically enables calculation of moments of the form E[X'^f.X]^], where m 
and n are at most equal to the number of orthogonal polynomials. This makes an important problem of 
constructing asymptotic confidence intervals for unknown parameters operational. 

Statistical test for Fisher-Snedecor distributional assumptions constructed in this paper is based on the 
generalized method of moments. This procedure is implied by the Stein equation (see Schoutens [64]) 
for the Fisher-Snedecor diffusion and on the results about the Pearson family of distributions of Bon- 
temps and Meddahi [21]. Beside these general results, Bontemps and Meddahi proved the robustness of 
Hermite polynomials for testing normality by using the Stein equation for the normal distribution (see 
[22]). However, we verified that Fisher-Snedecor polynomials are not robust test functions against the 
parameter uncertainty in our model (see Subsection 6.3.). 

Contents and main results. We focus here on statistical analysis of the ergodic diffusion process 
with Fisher-Snedecor invariant distributions. The paper is structured as follows. In Section 2 we present 
well known information on Fisher-Snedecor distribution and its moments. In Section 3 we define the 
Fisher-Snedecor diffusion and present its most important properties. Spectral representation of the tran- 
sition density of Fisher-Snedecor diffusion is developed in Section 4. In Section 5 we present method of 
moments estimators for parameters of this process and establish their asymptotic properties: consistency 
in Theorem 15.21 fi) and asymptotic normality in Theorem 15 . 21 f ii)- f iv V It is important to point out that 
orthonormality of Fisher-Snedecor polynomials and the closed form expression for spectral representation 
of transition density made possible calculation of explicit form of the limiting covariance matrix of the 
bivariate estimator of parameters of invariant distribution (see Theorem [521 (ii)). Section 6 develops 
methodology for testing statistical hypothesis about Fisher-Snedecor invariant distribution. The estab- 
lished method is based on the Stein equation for Fisher-Snedecor diffusion and the moment condition 
closely related to the Fisher-Snedecor polynomials. In Appendix A we present properties of the finite 
orthogonal system of Fisher-Snedecor polynomials. Appendix B contains well known facts about Gauss 
hypergeometric functions, while in Appendix C we present classification of boundaries of the state space 
of Fisher-Snedecor diffusion. 



2 General information about Fisher-Snedecor distribution 

Random variable X has Fisher-Snedecor distribution with a > and /3 > degrees of freedom, i.e. 
X ~ TS{a,l3), if its probabihty density function is given by 

where B{-,-) is the standard Beta function. 
Remark 2.1. While a more general density 

W2 (ax) 2 / N __ \aa;+«>'' ^ ^ ax+w ' ' 



might be of interest as well, in this paper we study the Fisher-Snedecor distribution with w = p. 

Our distribution belongs to the Pearson family of continuous distributions and is also known as Pearson 
type VI distribution (see Pearson [SI])- In particular, it follows that the tail of the Fisher-Snedecor 
distribution with density ([9]) decrease like a;~'-^+'^'^-', and so this distribution is heavy-tailed. 
Moment of the n-th order of J^S{a, (3) distribution is given by the expression 

''f''^-^nLr(/^-2^)-UJ r(|)r(f) ' ^"^ ' '^- 

Furthermore, moments could be calculated according to the recurrence relation 

In particular, expectation and variance are: 

Remark 2.2. Fisher-Snedecor distribution with positive integer degrees of freedom is frequently used 
in statistics (e.g. in analysis of variance). In particular, if x™ and Xm ^'^'^ independent chi-square 
random variables with n e N and ?7i e N degrees of freedom, respectively, then the random variable 
[(Xra/"-)/(Xm/"^)] has Fisher-Snedecor distribution FS{n,m). 



3 Fisher-Snedecor diffusion process 

Fisher-Snedecor diffusion is a solution of the non-linear stochastic differential equation 



dXt = -e (Xt - -^\ dt + J—^—Xt{aXt + /3) dWu t > 0, (10) 

(which is equation ([7]) for /i = /3/{/3 — 2)). The infinitesimal parameters, i.e. the drift parameter ^{x) 
and the diffusion parameter (t(x), of the Fisher-Snedecor diffusion (|10p are respectively given by 



^,{x) ^ -e { X ~ ^y a{x)^J-^^^x{ax + (3), /3 > 2, (11) 



where the restriction /3 > 2 ensures simultaneously that fj. = /3 / {f3 — 2) is positive and that 6 > 0, so that 

/i is the stationary mean. 

For a; > the corresponding scale density is 

s{x) = exp f- r ^^ du\ = x-^ {ax + P)^+^-\ (12) 

while the speed density is 

"^(^) = ^Tl^TT = ^^^-^-'{c.x + p)-^-^. (13) 

a'^[x)s{x) 20 

For any positive a and /3 the speed density (1131) is integrable on the diffusion state space, i.e. 

oo 

/m(.)ci. = ^fc^a'tr^i?(|,f)=A/<oo. (14) 



However, for the scale density s{x) we have 

s{x) dx = oo, Vxo € (0, oo), 



which ensures that starting from the arbitrary point xq from the interior of the diffusion state space the 
boundary oo almost surely cannot be attained (cf. Ai't-Sahalia [3]). A similar statement holds for the 
boundary 0, i.e. provided that diffusion starts from the arbitrary point xq G (0, oo), the boundary 
almost surely cannot be attained if and only if 

J.ix)dx = oo ^ a>2. 



We restrict ourselves to this case a > 2 (but note that in the opposite case a stationary diffusion may 
also be constructed, subject to instantaneous reflection at 0). 

For any positive a and /3 the stochastic differential equation (jlOp admits a unique strong Markovian 
solution {Xi, < i < T} with time- homogenous transition densities, since it satisfies the following 
sufficient conditions given by Ai't-Sahalia @]: 

(CI) the drift coefficient /i(a;) and the diffusion coefficient a{x) given by expressions (|11|) arc continuously 
differentiable in x on (0,oo) and cr'^{x) is strictly positive for all x G (0, oo), 

(C2) the speed measure p^ has the property given by ([Til) . 

According to Ai't-Sahalia |3|, these conditions are considerably less restrictive than the global Lipschitz 
and the linear growth conditions which are usually imposed on drift and diffusion coefficients to obtain 
existence and uniqueness of a strong solution (see Mikosch [60.)- 

For any a > 2 and /3 > 2 the Fisher-Snedecor diffusion is ergodic (see for example Genon-Catalot et 
al. [33], or S0rensen [7D]). If furthermore Xq ^ J^S{a,(3), then the Fisher-Snedecor diffusion is strictly 
stationary. For /3 > 2 the conditional expectation satisfy 

E [X,+t\Xs - x] - xe-«* + ^(1 - e-^*), 

and if /3 > 4, i.e. if the invariant distribution has finite variance, the autocorrelation function is given by 

p{t) = Coir iX,+t,Xs)=e-'>\ i > 0, s>0 (15) 

(see Bibby et al. [H], Theorem 2.3. (iii)). 



Remark 3.1. According to the general result by Genon-Catalot et al. (see [33], Corollary 2.1), the 
Fisher- Snedecor diffusion is an a-mixing process with an exponentially decaying rate, i.e. 

ax{t)= sup aiJ",, 7""+^) <]e-^\ 6 > 0, 

where 

a(J'^,J'"+*) = sup \P{Ar\B)-P{A)P{B)\, 

AeTs = cr{Xu, u<s}, Be J"''+* = cr{Xu, u>s + t}. 
Indeed, Fisher-Snedecor diffusion satisfies the sufficient properties by Genon-Catalot et al.: 

(i) the drift coefficient /i(x) and the squared diffusion coefficient (j'^{x), given by expressions (fTTj) . 
satisfy condition (CI). Furthermore, there exists a strictly positive constant 

a{l3-2) 

such that 

\^i{x)\ < K {I + \x\) and a^ix) < K {I + x^). 

(ii) The speed density m(a;) given by (fT3|) is integrable on the diffusion state space (see expression HH)). 
For a > 2 the scale density s(x) given by ((T^) is non- integrable in the neighborhood of boundary 
points and oo. 

(iii) The random variable Xq has (Fisher-Snedecor) density function which is proportional to the speed 

density m{x), i.e. 

, , , m(x) , , 

H^)^ ^^I(o,oo)(a;)- 

(iv) The product of the diffusion coefficient a{x^ and the speed density m(a;) converge to as x — > 
and a; — > oo, i.e. 

lim (T(2;)m(a;) — lim cr(a;)m(x) = 0. 

a:— ^0 a:— foo 

(v) If we define 



,, ,, , 2[i(x) /3 [a(x - 1) -I- 1] ex{ax + l3) ,^ , 

^{x)=a'ix)-^^ = ^ -L—^J^-—^^ x€0,oo, 

a{x) x[ax + p) y a(p — 2) 



then it follows 



1 a(/?-2) ,. ^xjax + (3) 

x^o j{x) y 13^9 x^o a{x - f ) + f 



J_^ / a(/?-2) ^xjax + (3) _ 1 1/3-2 ^ 

J^ ^^x) " Y /320 .^iSo a(a; - 1) + 1 " /? V 9 °^- 

Detailed exposition of the theory of mixing processes and additional information on corresponding central 
limit theorems are given by Doukhan |30j . 



Remark 3.2. The autonomous stochastic differential equation pll|) is a special case of the generalized 
non-linear mean reversion Ai't-Sahalia model 



dXt = (c-iX^' + CO + ciXt + C2Xf) dt + y/fjg + ajXt + cr|Xf dWt 
for following parameter values: 

c-i =0, erg = 0, 



/3-2 ^ "' "^1 — a(/3-2) 

-0 < 0, aj - - 

C2 = 0, 5 = 2. 



ci = -61 < 0, a| = ^ > 0, 



The generalized non-linear mean reversion Ai't-Sahalia model is frequently used for interest rates modeling. 
For more details about this model see [3] or |39) . 

Remark 3.3. The Lamperti transform {Yt, i > 0} of the Fisher-Snedecor diffusion is a solution of the 
stochastic differential equation 



dY,^^^^±^^^^t.nl.[Y,^Ji-^yt + dW,, t>0 

with unit diffusion parameter. According to lacus |39j . Lamperti transforms of diffusion processes are 
recommended for use in simulation studies. 



4 Spectral representation of the transition density of Fisher- 
Snedecor diffusion 

4.1 Fokker-Planck equation 

In this section we present the spectral representation of the transition density 

p = pix,t) =pix;xQ,t)^ —P{Xt<x\XQ^xo)., x>0, t>0, (16) 

dx 

of the Fisher-Snedecor diffusion in terms of solutions of the corresponding Sturm-Liouville equation. Note 
that the corresponding Sturm-Liouville operator is closely related to the infinitesimal generator of the 
Fisher-Snedecor diffusion. 

According to Karlin and Taylor (42] , the transition density ([T6|) is the principal solution of the Fokker- 
Planck or the forward Kolmogorov equation 

and its Laplace transform is known explicitly (see (p9)) ). 

Inverting the Laplace transform of the transition density yields its spectral decomposition - see for example 

the classical books by Titchmarsh [75], Karlin and Taylor [42 and Ito and McKean [40], as well as the 

paper by McKean ^S] (who call this " eigendifferential expansion"). 

The discrete part of the spectrum can be treated in a unified manner for all Pearson diffusions - see Wong 

[74] and Forman and Sorensen [32j . The exact formulas for the continuous part of the spectrum require 



however computations which need to be performed on a case by case basis. Some of such computations 
for special cases of Pearson diffusions were performed recently, motivated by applications in finance, by 
Linetsky [551 [531 [S3] , Davydov and Linetsky j^H] and Carr et al. [59^. We treated here the Fisher-Snedecor 
case, and we chose, for self containedness, to include a detailed account of the inversion of the Laplace 
transform of the corresponding transition density. 



4.2 Sturm-Liouville equation 

The infinitesimal generator of the Fisher-Snedecor diffusion is given by the expression 

/3 



a(/3-2) 

26* 

a(/3-2) 



(17) 



x{ax + f)f{x)-e\x 



P 



fix), x>0, 



acting on the domain 

D{g) = {fe L2((0,oo),fs(x)) nC2((0,oo)) : Gf e L\{0,^),fs{x)), 



lim 



fix) 



limlM^O 



>o s{x) K-s-oo s(a;) 

where s{x) is the scale density p^ . The negative of the infinitesimal generator G is called the Sturm- 
Liouville operator. The associated Sturm-Liouville differential equation {—Gf){x) — Xf{x) for the Fisher- 
Snedecor diffusion takes the form 



^^^^ x{ax + I3)f"{x) -9(^x- j^^ fix) + Xf{x) = 0, 



(18) 



where A > is the spectral parameter (see Karlin and Taylor [42j and Linetsky 
Remark 4.1. We consider also the equation 

Gfix) = sf{x), s > 0, 



(19) 



(with s — —A, where A is the spectral parameter from the equation I\l8\i ]. According to Ito and McKean 
(see [ID], page 128) or Borodin and Salminen (see [23], page 18), for any s > it admits two positive, 
linearly independent and increasing/decreasing solutions with Wronskian proportional to the scale density, 
called fundamental solutions. We will denote these solutions respectively by fi{x) ~ fi{x, s) and f4{x) = 
fi{x,s) (see (|24p and ([26| ). They intervene in the expression of Green's function, as well as in several 
first-passage problems. 




Figure 1. Graphs of fundamental solutions of equation l|19|l for o = 5, /3 = 20 and 9 = 0.05. 



10 



The Sturm-Liouville equation p5)) may be transformed by using the substitution 

/(x) ^v{y{x)) =v(-x\, 
to an equation of the hypergeometric type (see Appendix B): 

y{y + iy'{y) + (^| + (^1 - y^ v' {y) + X*v{y) ^ 0, (20) 



~ 261 • 

Denote the roots of the quadratic equation z^ + ^z + A* = by 



where A — — je 



.. = ... = -^±A„ A.A. = J^-M^. (21) 

± ±,A ^ A, ^ V 16 261 ^ ' 

According to Titchmarsh (see [72], Example 4.18., page 100), two linearly independent solutions of 
equation (^0]) for \y\ < 1 are 

V2{y) = y^~^2Fi(u+,u^;2~^:-yy a > 2, a ^ {2(to + 1), m G N}, 
where u± = u±.a = 1 — % + z±. Furthermore, two linearly independent solutions for \y\ > 1 are 

W3,4(y) = ?/"^^ 2i^i ( %,Wt;1t2Aa; — 

where the upper sign refers to the solution v^^y) and the lower sign refers to the solution V4{y). 

Remark 4.2. Cf. Abramowitz and Stegun (see [1], Section 15, Expression 15.5.21.), the existence of the 
solution V2{y) follows from the condition a ^ {2(m + 1), m G N} imposed on the value of parameter a. 
However, if a > 2 is an even integer, then the corresponding solution is given by much more complicated 
expression in terms of the digamma function ipl-) which is hardly evaluated in explicit calculations. 

Remark 4.3. According to Luke (see ^5, Section 3.9., Expression (1); see also Slater 68j, Section 1.8.1., 
Expression 1.8.1.11.), the solution vi{y) valid in the region \y\ < 1 can be analytically continued to the 
whole complex plane cut along the interval (— cx),— 1] (see Appendix B). This continuation is provided 
by the following expression which represents the solution vi {y) as the linear combination of the solutions 
V3{x) and V4{x), i.e. 

vi{x) = BxV3{x) + Axviix) (22) 

where 

_ r(f)r(2A,) _ r(f)r(-2A,) 

^"^'"r(z+,.)r(i -._,.)' ^ - "^^ - r (._,.) r(i -«+,.)• ^^^> 

Similarly, solutions v^ly) and ^4(2/), valid in the region \y\ > 1, could be analytically continued to the 
whole complex plane cut along the negative real axis. For details on the corresponding continuation 
formulas we refer to Luke (see [55], Section 3.9., Expressions (3) and (4)). This discussion shows that 
we are dealing with functions which are all analytic on the complex plane cutted along the subset of the 
negative real axis. 

Changing back variables, the solutions of Sturm-Liouville equation (IT5)) are 

/i(x)=/i(a;,-A) ^ ,F,(^z+,z_;^;-^xy (24) 
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/2(x)=/2(x,-A)^(^a;) \fi(^u+,u_;2-|;-^x), (25) 



a 



> 2, a(^{2{m + l),meN}, 



f3A{x)^f3A{x,~X)^(^x] ^ 2Fi(z^,u^;lT2Ay,—^Y (26) 

where the upper sign refers to the solution fsix, —A) and the lower sign refers to the solution fi{x, —A). 
Due to the analytic continuation of the hypergeometric functions, solutions fi{x,~X), /3(x, — A) and 
/4(x, — A) are all analytic on the state space of the Fisher-Snedecor diffusion. Therefore, fundamental 
solutions of equation (|19p on (0, cx)) could be identified with the increasing solution /i(x, s) and the 
decreasing solution fi{x,s), see Figure 1. 

Remark 4.4. According to Borodin and Salminen ([5F, page 19) and Buchholz |5T, Wronskian of 
linearly independent solutions fi{x,s) and f4{x,s) of Sturni-Liouville equation (1181) is proportional to 
the scale density ([T2|), i.e. 

WiUJi) = W,{UJi){x) - WsiUJi)ixo)six), (27) 

where the value a;o > could be chosen so that fi(xo,s), /{(xo,s), f4{xo,s) and /^(x^s) or the limit 
lim:r_j.^„ W {f4{x, s), fi{x, s)) are easy to calculate (see Titchmarsh [72]). According to expression (|C.3p 
from Appendix C, Wronskian of solutions fi{x, s) and /i(a;, s) is 



Wsih, h) = 2s{x)Bsa'-^r^ ^g + ^^^^ = 2B{x)a^~^r^B,As, (28) 

where Bg is given by (|23)) and As by (|2T|) for s = —A. 



4.3 Spectral representation of transition density 

As well-known (see e.g. Borodin and Salminen [53], page 19), the Laplace transform of transition density 
p(x;xo,t), also known as the resolvent kernel or Green's function 

oo 

Gs{xq,x) = / e'''*p{x;xo,t)dt, 



admits the explicit representation 

vnix] 

Gs(a;o,x) = — —- 'ipsixo A x)ipsixo\/ x), (29) 

where m(-) is the speed density, '0s(-) and ips{-) are the fundamental solutions of the equation (IT^ . and 
■Wsi'PsTi's) is their Wronskian (cf. Linetsky [52]) with respect to the scale density s{x), i.e. 

Ws{^s,'ips) = --^ {ip',ix)ips{x) - ipsix)(p'g{x)) = ' /[ " ■ 
S[X) s[x) 

Therefore, if it is possible to determine the explicit form of the fundamental solutions, then the transition 
density can be obtained by the Laplace inversion formula 

c+ioo 

p{x;xo,t) = — / e^*Gs(a;o,x)ds, i > 0, (30) 
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where s > is the spectral parameter from equation (J19p and where integration is performed afong a line 
Re(s) — c, c > 0, leaving all the singularities of the Green's function to its left. This approach implies the 
importance of understanding the nature of the spectrum of the Sturm-Liouville operator {—G) . Therefore, 
clarification of spectral properties regarding operator {—G) for Fisher-Snedecor diffusion is given in the 
following remark. 



Remark 4.5. According to the Theorem 16. II from Appendix C, boundaries of the state space of Fisher- 
Snedecor diffusion are classified as follows: is non-oscillatory boundary (regular for a < 2 and entrance 
for a > 2), while oo is oscillatory /non-oscillatory singular boundary (natural for all a > 0) with unique 
cutoff 

^=8(F^' ^>'- ^''^ 

Furthermore, oo is non-oscillatory for A = A (see Theorem I6.1l fii)). According to this classification 
of boundaries of the diffusion state space, the Sturm-Liouville operator {—G) has a finite set of simple 
eigenvalues 

in [0, A]. Therefore, the discrete spectrum of the operator (—G) is the finite set (Td{—G) = {A„, n = 0, . . . , [/3/4J }. 
Furthermore, the operator (—G) has the essential spectrum cressi~G) = [A, oo) which contains the purely 
absolutely continuous spectrum of multiplicity one in (A,oo), i.e. (Tad—G) = (A, oo). The elements of 
the absolutely continuous spectrum could be parameterized as 



'^^'^ " «rr2"''^~^"^' "'^'3'- 



/? > 2, (32) 



^ = ^ + l^^ = ^^[T:; + k'], /3>2, fc>0. (33) 



29P _ 20 //32 
13-2 ~ P-2 iy 16 

According to Linetsky's [55] spectral classification of one-dimensional diffusions based on the nature of the 
spectrum of Sturm-Liouville operator {—G), Fisher-Snedecor diffusion belongs to the spectral category 
IL 

According to Ito and McKean [JO] (see also Linetsky [55]), the general form of the spectral representation 
of transition density p^ for diffusions belonging to spectral category II is of the following form: 



N 



oo 



p{x]XQ,t)^m{x)^e ^"Vn(a;o)Vn(2;) +m(a;) /e ^Vl^^o, -A)(/3(x, -A) dA, (34) 

where A„, n e {0,...,A^}, are elements of the discrete spectrum (Jd{~G) and ^Pn{') are corresponding 
eigenfunctions normalized with respect to the speed density m{x), while A > A are elements of the 
absolutely continuous spectrum Uaci^G) and Lp{-^~X) are solutions of the equation (|18l) also normalized 
with respect to the speed density m(a;). 

Now we proceed to compute the spectral representation of the transition density for our process of interest 
following the general procedure from Ito and McKean [30] and Linetsky [52] and using the results of the 
analysis of the spectrum of the corresponding Sturm-Liouville operator (see Appendix C and Remark 
131). 
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Theorem 4.1. Spectral representation of the transition density of ergodic stationary diffusion with 
marginal Fisher- Snedecor distribution with parameters a > 2, a ^ {2(m + 1), to G N}, and (5 > 2 is 
of the form 

p{x]Xo,t) ^pd{x;xo,t) +pc{x;xo,t). (35) 

The discrete part of the spectral representation 



Pd 



{x;xa,t):^fs{x) ^^e^^"^ Fnixo) Fn{a 



(36) 



ra=0 



is given in terms of the eigenvalues A„ given by (1321) and the normalized Fisher- Snedecor polynomials 
Fn{-) given by (|A.3p . The continuous part of the spectral representation 



Pc{x;xQ,t) = fs(x) - 

TT 



'At 



fc(A)x 



(37) 



Bi 



Kt'l)r(-| + ^fc(A))r(f + | + zfc(A)) 



r(f) ^(l + 2^fc(A)) 



/i(xo,-A)/i(a;,-A)dA 



is given in terms of the elements A of the absolutely continuous spectrum of the operator {~G) given by 
P3|) . solution /i(-,— A) of the Sturm- Liouville equation ()18p given by (|24p and parameter k{X) = —iA\, 
where Ax is given in (j2ip . 

Proof. Since Wronskian of fundamental solutions fi{x, s) and f^ix, s) of equation ()19|) is given by expres- 
sion ([28l) . it is possible to obtain the explicit form of the Green's function for Fisher-Snedecor diffusion. 

In particular, 

m(a;) 



Gs{xo,x) 



2ai 



-Ib,a, 



■/i(a; Axo,s)/4(a; Va;o,s), 



(38) 



where m(-) is the speed density given by (J13p . /i(-,s) and /4(-,s) are non-normalized (with respect to 
the speed density) linearly independent solutions of the differential equation p9|) and 



Ws{fi,fi) = 2Bsa'-^r^\ TT; 



4. /3' , 5(/?-2) 



= 2a^-^r^BsAs 



16 29 

is their Wronskian with respect to the scale density ([T2|) . where Bg is given by (|23)) and A^ is given by 
(PT|) . Observe now the Green's function ([55| as the function of the complex variable s. We conclude that 
the term 



/3 , . /3^ 



^(^--) = M-i + Vi6 



20 



from i3s in expression p8|) has simple poles in 



— ~A„ — 



13-2 



n{l3-2n), n = 0, 



since for that value of s we have F (z+_a„) = T (— n), n = 0, . . . , [/3/4J. Note that these simple poles 
coincide with the negative simple eigenvalues of the Sturm-Liouville operator {—G) and that, since a > 2, 
a ^ {2(to + 1), to e N}, and (3 > 2, these are the only poles of the Green's function 
According to Abramowitz and Stegun [T], f4{x, s) = Csfi{x, s) + Dsf2{x, s), where 



C, 



F(l-f) F(1 + 2A,) 
F(u+,,)F(l-z_,,) 



and £>, 



F(f-1)F(1 + 2A,) 

F(z+,,)F(1-M_,,) 
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Therefore, regarding the fact that Green's functfon has simple poles at s = — A„, expression 
written in terms of solutions /i(-, s) and /2(-, s), i.e. 



Gs=-\„{xo,x) = m(x) 



r(-n)r(f + f 



ra-l 



a-fr^r(f)r(i + f 
r(f-i) 



(-1)" n(f + ^)/i(^0'^)/i(^'^)- 



can be 



(39) 



j=o 



+ - 



a 



-In 



/i(a;Aa;o,s)/2(a;Va;o,s) 



Due to supposed values of parameter a, the second term in expression p9p has no poles and therefore 
does not contribute to the residues of the Green's function. However, due to simple poles at s = — A^ 
and according to Titchmarsh (see [72], Chapter IV, Example 4.19.), the residues of Green's function at 
these poles are given by 



ReSs=-\^Gsixo,x) = m{x) 



(-1) 



2ri 



£_2nlr(f + f-n 



r/3-fr 



(f)r( 



l + f-nj ,=0 



n— 1 

n(|+.) 



X 2^1 



/3 a a 

-n, n ; — ; xn 

2' 2' 13 ' 



2^1 



~n, n - 



a a 



(40) 



2' 2' /3 

In expression (j40p we can recognize non-normalized (with respect to the speed density) Fisher-Snedecor 
polynomials given by formula (IA.3|) . 



Furthermore, Green's function 



has a branch point at s = — A = — 



ef3^ 



8(13-2) 



since for s < —A the 



gamma function r(z_|_ g) has the argument with non-zero imaginary part. Note that this branch point 
coincides with the negative cutoff between the discrete and the continuous part of the spectrum of the 
Sturm-Liouville operator {—G). The branch cut of discontinuity is placed from —A to — cx) on the negative 
part of the real axis and is parameterized as 

'/32 



-A 



29 



k' 



fc>0. 



13-2 V16 

After this analysis we can start with the evaluation of the inverse Laplace transform of Green's function 
(PS)) using the inversion formula (1501) . This procedure results in the spectral representation of transition 
density of Fisher-Snedecor diffusion. First we observe the Bromwich contour C as shown in the following 
figure. 
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Figure 2. The Bromwich contour C, simple poles, branch point and branch cut of Green's function. 
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Since simple poles (— A„) of the Green's function are placed inside the contour C, according to the Cauchy 
Residue Theorem it follows that 



2 



1 f LfJ 

— e'' Gs{xo,x)ds = y ReBs^-xy'Gs{xo,x). (41) 



c "=° 



On the other hand, the integral around the contour C is equal to the sum of the integral along the line 
AB, the integrals along the arcs BCD and HI A, the integrals along the lines DE and GH on each side 
of the branch cut, and the integral along the arc EFG around the branch point s = —A. 
Asymptotic behavior of the Green's function (I55|) as the radius r of the arc EFG around the branch point 
tends to zero implies that the integral along this arc vanishes. To verify this, first substitute s = re^'' — A, 
7 e (— TT, tt). Then the integral along the arc EFG has the following form: 

7r 

e'"Gs{xo,x)ds^~ie-^' / e''"'"'*G^e--,_A(a;o,a;)re*T ^7, (42) 

EFG — TT 

where 

Gre^._^{xo,x) = _, ^ ' h{xo A x,re'^ - A)f4{xo V x,re'^ - A). (43) 

2a^ 2/3 2Brg.^_^Arei^„A 

The expression under the integral vanishes as radius r — > if 

lini Greii~Aixo, x) < CX), 

which can be verified by observing the power series expansion of Greiy-A{xoix) around the point zero: 

Grei-i-A{xo,x) = c^^\ , „w „^+c^'^\r A^r x, ^ e' 2 ^ + c'^l , ^.. „ ^ e*'^r + 0(r^), (44) 

re '—l\\ u; / a,p,{xoAx),{xo\/x) a,!3,{xoAx),[xoVx) ^ aJ3,{xoAx),{xo\/x) ' V yi v y 

where quantities c^ a ^ ax) (x vxY * ^ {1:2,3}, are independent of r. Expression (|44)l implies that 
lim G^f.i-,_fXxQ-,x^ < oo. Therefore, the integral P^ vanishes as r — ^ 0. 



Furthermore, asymptotic behavior of the Green's function (|38)) as the radius R of the arcs BCD and 
HI A tends to infinity implies that the integrals along these arcs also vanish. Since these integrals are 
treated similarly, we explain the procedure regarding the arc BCD. Substitution s — Re^^ — A, where 
7 € (7r/2,7r), transforms the integral along the arc BCD in the following form: 

TT 

e'*Gs{xo, x) ds = ie-^* j Re'^''"'*GR^^-,_A{xo, x)e''' d-/, (45) 

BCD f 

where 



GRei-i-A{xo,x) ^m{x) 



n-f + V^^^f^e^* r f + H + WiH£^eH 



a-fr^r(t)r(^i + 2^M|_^e^* 

X /i(xoAa;,i?e*''-A)/4(a;oVa;,i?e^'^-A). (46) 

The expression under the integral on the right side of the expression (|^5t vanishes as the radius i? — > oo 
if 

lim Re^''''''^GRe^-,_t,{xo,x) = 0. 
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According to the formulas for asymptotic behavior of the Gauss hypergeometric function 2Fi{a,b;c; z) 
as two or more of its parameters tend to infinity (cf. Erdely [HI], page 77, expressions (16) and (17)) it 
fohows that 

fi{xo A X, Re'^ - A)/4(a;o V x, Re''' - A) = 



i^flLn 



» + |3-2 



» + )3-2 



r)£ + l I — ^ a,p,(xo/\x)' V ' a,p,(xo/\x)' \ a,p,{xoyx)> \ a,ti,{xayx)' 



ru + f + J^^^e^Mr i + V^^^^^^ 



r(t + | + Ji^i|^e^i r i-f-l + v/^^^e'Mr f-l + J^^^e^^ 



26* 



-a/^ 



i?(^ - 2) "^ ^ %,/3,(xoVa;) 



a,/3,(a:oA2;) 



V 



fiO-2) 



,,3,(xoA2;) 



1 + 



i?e'T 



(47) 



for large values of i?, where 

2a /4a 

fca,/3,(a:oAa;) = 1 + -^ (xq A x) + J ^ (^^O A x) (1 + (xq A x)) > 1, 

ka,0.Xxoyx) = 1 + ^(a;o Vx) + W— (xo Va;)(l + {xoV x)) > 1 
and fc„^^ (2;(,/^2.) < fca,^,(a:oV2;)- From expression (|47|) it follows that 

i?e^'='^*Gfle.T-A(xo,x) = 



^ ^(1-fc-i .O^ll + 'fc'^r A ^)'''^(l-fc~«r vO^(l + fc"«^ w O^'^^x 

„£ I 1 , — V a.pAxoAx)' y a,p.(xQ/\x)' \ q,/3.{xo Va;)' v Q,p,(a;oVa:) ' 

Z 2 ^-^ , /TT 



261 



R{p - 2) 



xi?e 



e '2 



fle'^'t 



r(i + f 



J?.(;3-2) g 
~^5 — ^ ^ 



Ji / fl(/3-2) .^ 
4 V 29 '^ 



r ( 1 - H _ £ 
^ ' ^ 2 4 



20 "^ I ^ 2" 4 



fl(/g-2) „i^ 



20 



("^Q,/9,(a;oAa;)"^Q,^,(2;oVa:)) * 



*a,/3,(a;oAa;) 



V^ 



(£-2) , 



'(*) //fc, 



*Q,/3,(2;oAa:) 
*a,,9,(2;oV3:) 



*v^ 



■(i) 



+ e"^(fc„,0,(^„^^)A:„,/3,(^,v:r))-' V 28 ™H2;(fc^^^_(^^^^^/s^^^_(^^^^j) 'V 28 '^'"UJ 



^ + ^'i^ 



for large values of R. Since fcQ,^.(2;oAz:) < ka^fi^(xoyx) and since for 7r/2 < 7 < tt we know that cos (^) > 
and cos (7) < it follows that i?e^^ ^Gjif^i-,_f^{x{)Tx) tends to zero as i? -> 00. Therefore, the integral 
(|15)l vanishes as i? — ?► 00. The same procedure is used for verifying that the integral along the arc HI A 
also vanishes as i? — ?> 00. Problem observed in a recent paper by Carr, Linetsky and Mendoza [59_ is 
treated with the similar technique. 

Finally, we conclude that the inversion formula pop for the Green's function (pS)) reduces to the sum of 
the integral along the line Re(a:;) = c, c > 0, and the integral of the jump across the branch cut, i.e. 



c+ioo 



2-Ki 



e'*Gs{xo,x) 



2ni 



2TTi 



e''Gs{xo,x)ds + — / e^* (G,(/i,/4)-G,(/i,/4)) ds, (48) 
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where Gs{fi,f4) is the complex conjugate of the Green's function ([55)) . From expressions (|1T|) and 
it follows: 

c+ioo I |. I -A 

^ / e^*G,(xo,a:)ds = ^Res3_A„e^'G,(a;o,a;) + ^ / e^* (Gs{h, h) - Gs{h, h)) ds. (49) 



n=0 



Now it remains to determine the explicit form of the expression Gg (/i , /4) — Gg (/i , fi) for the jump across 
the branch cut. Since /i(-,s) = /i(-,s), fi[-,s) = fz[-,s) and /i(-,s) = Bfj,{-,s) + AJi{-,s) (see Remark 
1221 ). for s = — A the jump across the branch cut is given by expression 





G-x{xa,x) - G'_A(a;o,a;) = 


m{x)2ik{\) 


a°4 V' r (-f + ik{\)) r (f + 1 + zfc(A)) 


r(f)r(i + 2ifc(A)) 



(50) 



/i(a;o)/i(2^)- 



Substitution of the expressions (PH)) and (ISUl) in the relation (P^)) results in the spectral representation of 
the transition density of Fisher-Snedecor diffusion. However, Fisher-Snedecor polynomials and solutions 
of the Sturm-Liouville equation (fT8| for A > A are not normalized with respect to the Fisher-Snedecor 
density ([9|). Now the uniqueness of solutions of the Sturm-Liouville equation up to a constant factor 
ensures that instead of the solution /i(x, — A) we can use the solution fi{x,—X) — J wz2 hi^^~^)- 
This procedure results in the expression for the spectral representation of transition density in terms of 
functions normalized with respect to the Fisher-Snedecor density, i.e. the result is exactly the expression 
p{x;xo,t) = pd{x;xo,t) + pc{x;xo,t), where pd{x;xo,t) is given by ([36]) and pc{x;xo,t) by ^7^. D 

Remark 4.6. In the previous theorem we supposed that a > 2, a ^ {2(m + 1), m E N}. If this is not 
the case, i.e. if a > 2 is an even integer, then in pop expression for the solution /2(a;, s) discussed in the 
Remark 14.21 should be used. However, due to its complicated form it is hard to expect that a closed 
form expression for spectral representation of transition density is obtainable in such an elegant form as 
in the Theorem 14. II 



Remark 4.7. The general approach to the spectral theory of one-dimensional diffusions is treated in 
a number of papers by Linetsky (see [S5J [531 [Ml El]). Beside the concise and informative introduction 
to basics of the spectral theory, he gives the procedure for calculation of the spectral representation of 
the transition density for one-dimensional diffusions (see [SH [Sl]). We already used his approach for 
transition densities of reciprocal gamma and Student diffusion (see [49l [50l [13]). Since these processes, 
together with the Fisher-Snedecor diffusion, make a class of heavy-tailed Pearson diffusions, it is natural 
that spectral representations of their transition densities have similar forms. Namely, for all three heavy- 
tailed Pearson diffusions spectral representations of transition densities are given in terms of the finite 
sum of residues and an integral of the jump across the branch cut of the corresponding Green's function. 
However, the differences between mentioned spectral representations are generated by different types of 
the corresponding Sturm-Liouville equations. In particular, Sturm-Liouville equation related to reciprocal 
gamma diffusion can be transformed to the Whittaker equation, while the Sturm-Liouville equations for 
Student and Fisher-Snedecor diffusions can be transformed to the hypergeometric equations. According to 
the qualitative nature of the spectrum of the corresponding Sturm-Liouville operator {—G), the continuous 
part of the spectral representation of transition density of reciprocal gamma diffusion is written in terms 
of Whittaker functions or hypergeometric function 2-Fo(ai b; ; •) (see |S31 ISll I3S1 [Zl] ) ; while the continuous 
part of the spectral representation for Student and Fisher-Snedecor diffusions is written in terms of the 
hypergeometric function 2^1 (a, b; c; •), see Appendix B. 
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Remark 4.8. Fisher-Snedecor polynomials Fn{x), n — 0,. . . , [/3/4J, and functions fi{x, —A) for A > A 
belong to orthogonal subspaces of the Hilbert space V. ~ L^((0,oo),fs(-)), i.e. 

Fr,ix)enpp, nG{0,l,...,L/3/4j}, /i(x,-A)eH„„ VA > A. 

Here Upp denotes the subspace of the Hilbert space L~((0, oo), fs(-)) containing functions having only 
the pure point spectral measure, while Hac denotes the subspace of the Hilbert space L^((0, oo), fs(-)) 
containing functions having only the spectral measure which is absolutely continuous with respect to the 
Lebesgue measure (cf. Reed and Simon [B^ or Linetsky [S^])- From here it follows that Fisher-Snedecor 
polynomials and functions /i(x, — A) for A > A are orthogonal with respect to the invariant density of 
Fisher-Snedecor diffusion, i.e. 

oo 

' Fn{x)fi{x,~\)fsix)dx = 0. (51) 

This orthogonality relation plays a crucial role in the statistical analysis of this diffusion process. 

Remark 4.9. Two-dimensional density of the stationary diffusion process X = {Xt, t > 0} that solves 
stochastic differential equation ([TU]) and has marginal density function is given by the expression 

P{x,y,t) = -^^P{Xs+t <x,Xs< y) = fs{y)p{x;y,t) = fs{y) {pd{x;y,t) + pc{x;y,t)) (52) 

where Pdix; y, t) is given by (P5)) while Pc{x] y, t) is given by ((57)1 . 



5 Parameter estimation for Fisher-Snedecor diffusion 

5.1 Estimation of autocorrelation parameter 9 

The autocorrelation parameter 9 > Q will be estimated under the assumption that parameters a and /3 
of the Fisher-Snedecor diffusion are known. The ergodicity of the Fisher-Snedecor diffusion is ensured by 
the assumption a > 2 while the existence of the first and the second moment of the invariant distribution 
is ensured by the assumption ji > A. 

Let us consider the sample Xi, . . . , Ar„ from the ergodic stationary Fisher-Snedecor diffusion and the cor- 
responding sample of paired observations {Xi,Xt+i), (X2, ^4+2), ■ • ■ , [Xn-t^ A'„), where t < n. Empirical 
counterpart of the autocorrelation function is given by the Pearson's sample correlation coefficient 

Pn{t) = ^ I - , (53) 

n-t l^i=l ^i \n-t l^i=l ^y V "-* ^i=l t+i \n~-t l^i=\ ^t+t J 

where the term in the numerator represents the empirical covariance of random variables Xs and Xt+s, 
while the term in the denominator represents the product of the empirical standard deviations of random 
variables Xg and Xt^g, respectively. 

Expression (|53p implies that there exists a continuous function g: M^ \ {{x,y, z,u,v) : x^ ~ z or y^ — 
u} — > M defined by the expression 

v-xy 
g[x,y,z,u,v} 



^ {z - x'^){u - y'^y 
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such that 

/ n—t n—t n—t n—t n—t 

\ z— 1 z— 1 -i— 1 -i— 1 i— 1 

Since Xi, . . . , X^ is the sample from the ergodic stationary Markov process, according to ergodic theorem 
for stationary sequences (see KarHn and Taylor [51], Theorems 5.4. and 5.6.), it follows that 

(^ n—t ^ 7i~t ^ n—t ^7i—t ^ n—t \ 

^ E ^- ^ E ^*+- ^ E ^'' ^ E ^'+" ^ E ^^^*+^ (54) 

n — f ■'^-^ n — i -^-^ n — t ^-^ n — t ^-^ n — t ^-^ ] 

1=1 i=l i=l 1=1 i=l / 

is the P-consistcnt estimator of E [(Xq, Xj, Xq , Xj'^, XgXf )] . 

FroiTL expression (J15p for the autocorrelation function p(i) of the Fisher-Snedecor diffusion we see that it 
takes only positive values. However, according to Hassani (cf. [37 , Theorem 2.1.) the sum of the sample 
autocorrelation function at lag i > 1 is always —1/2 for any stationary time series with arbitrary length 
71 > 2, i.e. 

n— 1 -. 

SacS = X!?"(*) "" ^2' 

Therefore, instead of the Pearson's sample correlation coefhcient ([55]) we observe here its absolute value 
|p„(i)|- Since |pn(i)| is the continuous transformation of the P-consistent estimator ([54l) . according to 
continuity mapping theorem it is the P-consistent estimator of the autocorrelation function (|15p for any 
fixed i > 0, i.e. 

|Pn(i)| — > p(0; "^ ~^ °^' foi' ^^y fixed t > 0. (55) 

Therefore, here we observe the equation |/0„(t)| — p{t\ where p(t) = e^^*. Solving this estimation 
equation in terms of the unknown parameter Q results in the estimator of the autocorrelation parameter 
6*, i.e. 

^=-ilog|p„(i)|. (56) 

The continuity of the logarithmic function and the relation (j55p imply that 

1 _ pi 

log|p„(^)| — > logp(t), n — ^ oo, for any fixed i > 0, 



which means that d is the P-consistent estimator of the autocorrelation parameter 0, i.e. 



(57) 



5.2 Estimation of parameters a and /3 

In this section we present method of moments estimators of unknown parameters a and /3 of invari- 
ant distribution of ergodic Fisher-Snedecor diffusion and analyze their asymptotic properties under the 
assumption that d is the known value of the autocorrelation parameter. Estimation of the unknown pa- 
rameter (a, /3, Q) by the method based on the martingale estimation functions (cf. Forman and S0rensen 
[32]) is briefly discussed in Section [Ol However, in contrast to the method which will be developed here, 
their method for calculation of limiting covariance matrix in asymptotic normality framework seems to 
be more complicated. 

A new method for calculation oi E\X'^j^^X^\ based on the orthogonality of the normalized Fisher-Snedecor 
polynomials with respect to the density ([9|) will be introduced. The existence of Fisher-Snedecor poly- 
nomials F\(x) and ^2(2;) needed for mentioned calculations is ensured by the assumption that /3 > 8. 
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Equating the first and the second theoretical moments with the corresponding empirical counterparts 

rni = -J2^t and ^2 = - ^1 ^? (5^) 

and solving the resulting system of equations in terms of the unknown parameters a and (3 results in the 
method of moments estimators for these parameters. They are given by following expressions: 

« = _ ,^ _\ _2 , (59) 

7712(2 — 7771 ) ^ "T-l 

d=^. (60) 

The statement of the following proposition follows from the general result for the spectral representation 
of transition density for one-dimensional diffusions. This special case for Fisher-Snedecor diffusion is 
treated here in details to demonstrate the use of a closed form expression for spectral representation of 
transition density (see Theorem 14.11 '). Analogue results for reciprocal gamma and Student diffusions are 
presented in [49l [50] . 



Proposition 5.1. Fix i and j such that i,j G {!,..., [/3/4J} a77d suppose that j3 > 2{i + j), a > 2, 
a ^ {2(777 + 1), 777 e N} and 9 > 0. If {Xs+t,Xs) is a two-dimensional random vector with the density 
function (|52p and F„{x), n G {1, . . . , [/3/4J}, are orthonormal Fisher-Snedecor polynomials given by 
Rodrigues formula (|A.5|) , then 

E[F,{X,+t)Fj{X,)]=e-^^'d,j, 

where Xj is the eigenvalue given by p2p and Sij the standard Kronecker symbol. 

Proof. The two-dimensional density function p{x, y, t) given by (j52[) could be written in the following 
form: 

LIJ 

p(x,y,i) = ^e-*^"*n(.T)<I>„(2/) + - / e-'''u\X)-^{x,\)-^[y,\)d\ (61) 



" = 8/32 



where 





v\\) - 


--k{\) 


^K^' 2) r (-1 + *^W) r (f + 1 + ^MA)) 




r(f) r(i + 27fc(A)) 


and where 




*«(•) = H-)M-), 






*( 


,A) ^ fs(.)/i(-) = fs(-)/i(-,-A). 



Here Fn{-)^ n € {1,..., [/3/4J}, are Fisher-Snedecor polynomials and yV^(A) /i(-, — A) are solutions of 

the Sturm-Liouville equation (1181) for A > A, both normalized with respect to the invariant density of the 

Fisher-Snedecor diffusion. 

Orthogonality of Fisher-Snedecor polynomials Fn{-), n E {1, . . . , [/3/4J}, and functions /i(-, —A), A > A, 

(see the Remark 14. 8p together with the representation (pT|) of two-dimensional density ([5^ imply the 

following: 



E[F,iX,+t) Fj{Xs)] ^ / F,{x)F,{y)pix,y,t)dxdy 
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F,ix)F^ix)f5ix)dx\ \ / F,(y)F„(y)fs(2/)dy 








F,ix)h{x,X)fsix)dx\ F,iy)My,X)fs{y)dy\ dX 



-^'*5,,=exp^-0j^-|t^.,„ 



where, according to the Reniark l4.8[ the term related to the continuous part of the spectral representation 
psp disappears. In order to apply the Fubini's theorem in the above computation we used the property 
that the solution /i(z, —A) = 2-Fi(a, 6; c; —z) of the hypergeometric equation can be written in the form 

/i(z, -A) = C^{-z)- + C2{-z)-' + 0{{-z)—^) + Oll-z)-"-!), 

where Ci and C2 are constants and (a — b) is not an integer. D 



Expressions (|59p and (|60p imply that the estimators a and j3 are continuous transformations of the 
empirical counterparts of the first and the second moment given by ()58|) . In particular, a — gi(jni,m2) 
and /3 — (72(^1,7712). Therefore, asymptotic properties of the estimators a and /3 are directly implied 
by asymptotic properties of the empirical counterpsrts mi and 7772, which will be presented now. Since 
Fisher-Snedecor diffusion has the a-mixing property with the exponentially decaying rate (see Remark 
13. ip . in further analysis we can use the following central limit theorem for a-mixing sequences (see Hall 
and Heyde [36] or Billingsley [2^; for multidimensional version sec Genon-Catalot et al. [34]). 

Theorem 5.1. Let < S < 00. Suppose that X ~ (Xn, n G N) is a stationary and a- mixing sequence 
of random variables such that i?[|Xo|] < 00 and let /i, ...,/„ : R ^ R, m E N, be Borcl functions such 
that 

• S[|/fe(Xo)|2+^]<^ifO<(5<oo, 

• \fk{Xo)\ < c < 00 a s = 00, 



• E(a/fe(x)(n))^^ < 00 

n=l 



Let §„ = ^(/i(Xfe), . . . , U{Xk)). Then 



k=l 

00 CX) 

4 = Cov{MXo)Jj{Xo)) + Y. GoviMXo), MXk)) + ^ Cov(/,(Xfe), /,(Xo)) 

fe=i fc=i 

exists yi,j G {1, . . . , 777}. If the matrix E — [dfA^mxm) is positive semidefinite, then 

^ — ^^A/(0, E), 77 -> cx). 

Remark 5.1. As an alternative, the central limit theorem can be formulated in terms of the /3-mixing 
property of stationary diffusion. Beside the approach of Genon-Catalot et al. [Sj, /3-mixing property of 
Fisher-Snedecor diffusion can be proved by the methodology proposed by Abourashchi and Veretennikov. 
In recent papers (see [21 13]) these authors proved that reciprocal gamma and Student diffusion models 
posses the exponential bounds for the /3-mixing coefficient both in time and in space. The latter result 
is particulary valuable since it implies that functional central limit theorem can be also formulated for 
non-stationary reciprocal gamma and Student diffusion models. 
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Remark 5.2. An alternative approach to the central limit theorem for stationary diffusions is based on 
the theory of associated random systems. This approach is developed by Bulinski and Shashkin (see pS] . 
Corollary 2.30. and Theorem 2.29.) and takes its origin in the paper by Barbato (see [14,, Lemma 8) 
where he proved that the stationary diffusions satisfy association conditions. 

The ergodic theorem for stationary sequences (see Karlin and Taylor [41 or the Lemma 20.3. from 
Billingsley [19]) imply that {rni,rn2) is the P-consistent estimator of E[{Xt,Xf)], while Theorem 15.11 
implies that {rni,rn2) is asymptotically normal, i.e. 



y^{mi - E[rni],m2 ~ E[m2]) 4 A/'(0,S) 



(62) 



where 



E[mi] 



P 



E[m2] 



/32(a + 2) 



/?-2' ^-"^^ a(/3-2)(/3-4)- 
To obtain the explicit form of the (2 x 2) covariance matrix S we must first calculate the expectations 
E[Xs+tXs\, E[Xs+tXf\ and E[Xl^iXl], where Xg and Xg+t are random variables from the Fisher- 
Snedecor diffusion. These moments can be calculated by using the known expectations of the products 
of Fisher-Snedecor transformations of random variables Xg and Xs+t- In particular, we observe random 
variables F„(Xs) and Fn{Xsj^t), where n € {1,2} and F„(-) is normalized Fisher-Snedecor polynomial 
(jA.Sp . According to the Proposition 15. 1[ expectations of their products are given by 

E[F^{Xs+t)Fi{Xs)]=e-^\ 
E[F^{X,+t)F2{X,)]^Q, 



E[F2{X,+t)F2{X,)]=e 



'29J^t 



On the other hand, the same expectations can be calculated regarding the explicit definitions of the 
Fisher-Snedecor polynomials Fi{x) and ^2(2;) (see Appendix A, section A. 2.) and written in terms of 
the unknown expectations E[Xs+tXs], ElXs+tX"^] and E[X'^_^_f.X'^]. In particular. 



E[Fi{Xs+t)Fi{X,)] = 



«(/3-4) 
2/32(a-f/3-2) 



[ip-2fE[X,+tX,]~l3^], 



E[FiiX,+t)F2iXs)] 



'(^-2)(/3-4)(/3-8), 3 



4a/33(a + /? - 2) V {a + 2)(a -f /? - 4) 



a^(/3-2)(/3-4)(/3-6)x 



xE[Xs+tX^] 



Aa/S^a + 2){a + 13 ~ 2) 



-9t 



a^p^{a + 2){p~6) 



v{f3 - 2){(3 - 8) 



S [F2(X,+,)^2(X.)] - ^^ ^ 2)(a + /3 ^ 2)(a -H /3 - 4) 



P-2 P-2 

'a\f3~Af{P~QfE[Xl^,Xl] 



Ka + 2f 



A){a + fi-2) 



-et 



/34(a + 2)2(/3-6) 
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a(/3-2)2 



(/3-2)2 



This approach results in the system of equations whose unknowns are expectations E [X^+f X^] , E [Xg+f X^ 
and E[X1^^X1]. Explicitly, these expectations are given by the following expressions: 



E [Xsj^t Xs\ 



E[X,+tX^] 



(/3-2)2 



2/32(q + /3-2) 
a(/3-4) 



-et 



P' 



/33(a + 2) 



a2(/3-2)2(/3-4)(;3-6) 



[4(a + /? - 2)6-"* + a(/3 - 6)] = E[xI+^X,], 
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pr^2 y2i _ 8(3\a + 2){a + /3 - 2)(a + /3 - 4) ,2, fel . 
L *+* *J a3(/3-2)(/3-4)2(/3-6)2(/3-8) 

, 8/34(a + 2)2(a + /3-2) _ej , /S^Ca + 2)^ 



a3(/3 - 2)2(/3 - 4)(^ - 6)2 a2(/3 - 2)2(/3 - 4)2 ' 

Elements cr?-, i,j € {1,2}, of the symnietric covariance matrix S now can be determined according to 
the formula for covariances given in the Theorem 15.11 In particular, 

2 2/32(a + /3-2) e^ + 1 



a(/?- 2)2(^-4) e^'-r 



4/32(a + 2)(a + /3^2) e^ + 1 _ 3 
a2(/3 - 2)2(/3 - 4)(/3 - 6) e» - 1 ~ ^^^i' 

/3^(a + 2) [(q + 4)(a + 6){J3 ~ 2)(/3 - 4) - a(a + 2)(/3 - 6)(;g - 8)] 
a3(/3-2)2(/3-4)2(/3-6)(/3-8) + 

16/3'*(a + 2)(a + /3-2)(a + ;3-4) 1 16^'*(a + 2)2(a + /3 - 2) 



a3(/3-2)(/3-4)2(/3-6)2(/3-8) g^^l^ - 1 ^^(/S - 2)2(/3 - 4)(/3 - 6)2 e« - 1 ' 

Since a > 2, a ^ {2(m +1), m G N}, /3 > 8 and 6* > 0, the covariance matrix S is positive definite. This 
concludes the analysis of asymptotic properties of the estimator (7711,7712). Now we can start with the 
analysis of asymptotic properties of the estimator ( a, /? J , given in the following theorem. 

Theorem 5.2. Let {Xt, t > 0} be the Fisher- Snedecor diffusion with the unknown parameter (a,/3), 
where a > 2, a <^ {2{m +1), m e N}, and /3 > 8. 

(i) (a, /3) — >■ [a, 13), i.e. {a,j3) is the P-consistent estimator of the unknown parameter {a, [3). 

(a) ^/n{a — a, /3 ~ 13) — > M{0, S(a, /3, 0)), where S(a, /3, 6) is the positive definite covariance matrix. 

r ^ -" 1 ~-'-/2 ^ ^ d 

(Hi) \/n S(a, /3,6') (S — a,/3 — /3) — > A/'(0, 1), where is the known value of the autocorrelation 

parameter. 

r ^ -1 -1/2 ^ '^ d -- 

(iv) ^/n 'S{a,/3,0) {a — a, (3 — (3) — > M{Q,1), where 9 is the P-consistent estimator ([5B)) 0/ the 

autocorrelation parameter 9. 

Proof. 

(i) Estimators mi and m2 given by (j58p are the P-consistent estimators of the first and the second 
moment, respectively. Since estimators a and /3 are their continuous transformations, i.e. 

(S,^) = (51 (^1,^2), 52(^1, "^2)), 

according to the Theorem 1.7.1. from Serfiing [5S], it follows that 

a^ gi{mi,rn2) ^ giiE[rni,m2]) ^ a, n -> 00, 



j3 = 52(^1,7712) -> 32 {E[mi,m2]) ^ P, n^ 00. 



From here it follows that 



^ p 

{a,P) ^ {a,l3), n-^oo, 



i.e. (S, (3) is the P-consistent estimator of the unknown parameter (a, f3). 
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(ii) Since a and /3 are the continuous transformations of mi and 7712, according to the muhivariate delta 
method (see Serfling [SS], Theorem 3. 3. A.), the estimator (a,/3) is asymptotically normal, i.e. 



y/7i{a - aj - 13) A Af {0, DSD^) 



dgi 



dxi 



where a = gi {E[mi, 7712]), P — g2 (£'[7^1, 7rt2]) and (2 x 2) matrix D — 
has the following form: 

a(a + 2)(/3-2)(3/3-8) a^il3 - 2){I3 - i) 



(63) 
,»e{i,2}, 



X— £^[rni ,m2] 



D 



2/3(/3 - 4) 
(/3-2)2 



2^2 




Explicit forms of the elements afja,(5,6) of the covariance matrix D'ED'^ = 'E{a,f3,9) are given 
by the following expressions: 



aUa,(3,0) = 



k(a + 2)(/3-2)(a + /3-2) 



2/3(^-4)3(^-6)2 
(a + 2)(-3072 + /3(6528 + /3(-4736 + /3(1548 + /3(-232 + 13/3))))) 



4/3(/3-4)3(a + /3-4) 
/3-8 



coth 



/?- 


-2 


0(/3- 


-4)\1 


/3- 


2 JJ 



coth 



(g + 2)(/3 - 2)(a + /3 - 2) [/3(/3 - 4)(3/3 - 16) - 32] 

crfn a, /3, 6* = ^ -77^ ^ coth , 

12^ ''^' ' 2(/3-4)2(/3-6) V2 

9. ^ .^ ;3^(/3-2)2(a + /3-2) ,/6'\ 

erf, a, /3, 6*) = — !^- P f- '- coth - , where coth , , „ . 

^^^ ''^' ' 2a(^-4 V2y' V2y e«-l 



ai,ia,p,e), 



(iii) Part (ii) of this theorem and the positive definiteness of the covariance matrix S(a, /3, 0) imply that 



[S(a, /3, 0)]-'^' V^ (a - a, /3 - /3) A AA(0, 1), 



(64) 



where I is the (2 x 2) identity matrix. Since (S, /3) is the P-consistent estimator of the parameter 

{a,/3), it follows that 

-, -1 /9 . „ 

(65) 



S(S,/3,0) ''^' [S(a,/3,0)] 



1/2 p 



If we multiply 



from the left by the matrix S(a, /3 

(see Serfling [65j, theorem 1.5.4.1.) imply that 

,P,e)Y'^' [S(a,/3, 

Since [E{a,l3,d)]^/^ [S(a,;9,6')]"^/^ 



[S(a,/3,0)]"'/' V^{a-aJ~l3) 
-1/2 



S(a,/3,6')]^/^ then dH]), dM]) and the Slutsky Lemma 



s(a,/3, 



11/2 



S(a,/3,0)]^/"[S(a,/3,0)] 



-1/2 



■(a-a,;3-;9)4AA(0,I). 



S(S,/3,0) 



= I, it follows that 
1/2 



(3-a,/3-;3) AAA(0,I). 
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(iv) Since 6 given by (|56p is the P-consistent estimator of the unknown autocorrelation parameter 9, 
the proof foUows by the same argument as the proof of (iii) . 

D 

Remark 5.3. Since [S(a,/?,6')]~ ' \/n ( a — a, ^ — /3) — > A/'(0,I), from expressions for elements of 

the asymptotic covariance matrix it follows that ( S, /3 ) is the y^-consistent estimator of the unknown 
parameter (a, j3). 

Remark 5.4. Estimators ([5^ and (|5U1) of parameters a and (3 are the method of moments estimators 
based on the empirical counterparts mi and m2 given by (j58p . It is interesting that if instead of m,2 the 
empirical counterpart of the second central moment 



1 " 

n ^ — ^ 



n 

is used, the method of moments results in the same estimator for parameter a, i.e. the result is the 
estimator (I59p . In particular, the use of M2 results in the estimator of parameter a which is given by 
the expression 

a = = ^ ■ (66) 

M2(2-TOi)-TO2(mi -1) 

Since M2 = m2 — Tn^ , substituting this in the expression (|66|) results in the expression (|59p . 

5.3 Martingale estimating functions 

In this section we briefly discuss estimation of the unknown parameter ip = (a, /3, 9) of the Fisher-Snedecor 
diffusion by the optimally weighted martingale estimating function approach due to Forman and S0rensen 
|32) . The result presented here is just a special case of the more general result for estimation of unknown 
parameters of Pearson diffusions, specialized for the Fisher-Snedecor diffusion. 
In particular, the martingale estimating function for the Fisher-Snedecor diffusion is given by 

n N 

where 

• Wj{-,ijj) are the optimal weight functions chosen to minimize the asymptotic variance of the esti- 
mator (cf. expressions (3.4) and (3.5) in [35]), 

• eigenvalues Xj given by expression (I32p depend on the unknown parameters /3 and 9, 

• Fisher-Snedecor polynomials Fj{-) defined by Rodrigues formula (IA.5I) depend on the unknown 
parameters a and /?. 

According to the Theorem 3.1. from [32], the system of martingale estimating equations Gnii/j) = 
provides consistent and asymptotically normal estimator ip — {a, (3,9) of the parameter ijj = {a, (3,9) 
of the ergodic Fisher-Snedecor diffusion with the restriction (3 > 2N on the value of the parameter (3, 
ensuring the existence of the first N Fisher-Snedecor polynomials, if the mapping -0 1— >■ ip — (9, fi,b2,bi,b()) 
is differentiable and the matrix df/dtp'^ has the full rank, i.e. the rank three. Here 

[(3 2 2(3 

^^i9,f^,b2,b„bo)^[9,j-^,j-^,-^^-^^,0 
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is the canonical parameter of the ergodic Fisher-Snedecor diffusion (cf. equation (2.1.)) and the corre- 
sponding matrix of partial derivatives 



dip 



n 2 2 4 Q 

'-' (/3-2)2 (^-2)2 a(0-2)2 " 

10 



has the rank three. Therefore, Theorem 3.1. from Forman and S0rensen [35] ensures the existence of 
consistent and asymptotically normal estimator of the parameter tfj = {a, /?, 9) of the Fisher-Snedecor 
diffusion. The procedure for calculating the explicit form of the asymptotic covariance matrix in the 
asymptotic normality framework is also given in the Theorem 3.1., but this procedure seems to be more 
complicated than the procedure based on the orthogonality of Fisher-Snedecor polynomials and related to 
method of moments estimators of the parameter (a, /3) of the invariant distribution of the Fisher-Snedecor 
diffusion (cf. Subsection 5.2.). 

In conclusion, martingale estimation function approach provides optimal estimators of the unknown 
parameters of the Fisher-Snedecor diffusion, in the sense that the variances of estimators are minimized, 
but the calculation of the asymptotic covariance matrix could be complicated. On the other hand, 
method of moments provides estimators in the explicit form and, since these estimators are continuous 
transformations of the empirical counterparts of the first and the second moment, there is an elegant 
procedure for calculation of the asymptotic covariance matrix in the asymptotic normality framework 
which makes this method more applicable and easier to use. 

6 Testing Fisher-Snedecor distributional assumptions 

6.1 Stein equation for Fisher-Snedecor diffusion 

The results of this section are based on the general theory developed by Barbour [15, and Schoutens [64] 
and derived under the following assumptions on values of parameters a and /3: 

a>2, a^ {2(m+l), me N}, ^ > 4. 

The general form of the Stein equation is 

h{x) - E[h{Z)] = Afh{x), (67) 

where h and fh are smooth functions, Z is the random variable with the known distribution and A is the 
Stein operator closely related to the distribution of the random variable Z. For every smooth function 
h the wanted solution of the equation (|571) is the function //, with the property that for each random 
variable Y the equality 

E[h{Y)] - E[h{Z)] = E[Afh{Y)] 

holds. According to Schoutens [64!, distribution of the random variable Y is close to the known distri- 
bution of the random variable Z if and only if E[Afh{Y)] is close to zero for a large class of functions 
fh and a proper choice of the Stein operator A. According to Barbour [15 , the best choice for the Stein 
operator A related to the invariant distribution of the diffusion process is the corresponding infinitesimal 
generator. 
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Substitution of the infinitesinial generator (I17|) of the Fisher-Snedecor diffusion and the random variable 
Xt ~ FS{a,(3) in the equation (|57)) results in the Stein-Markov equation 

20 / /3 



h{x) E[h{X,)] = -^^^ x{ax + P)g"{x) ~ 6 (^x ~ ^) 



5'(a:), 



where g(-) is the smooth function. Substitution of the function fh{x) — g'{x) into the equation 
results in the equation 

h{x) - E[h{Xt)] - 



29 



■x{ax + P)f',^{x)-e[x 



P 



h{x) 



(69) 



which will be called the Stein equation for Fisher-Snedecor diffusion. Considering Barbour's method, the 
corresponding Stein operator is given by the expression 



Ah{x) 



29 



a(/3-2) 



x{ax + P)f;^{x)-9[x 



(3 



p-2 



fh{x), X e (0,oo). 



(70) 



According to the general result from Schoutens [64^ for Pearson diffusions, solution f^ of the Stein equation 
([S9)) is of the form 



h{x) 



a'^{x)^B{x) 

~2 

a'^{x)^5{x) 



{h{y)~E[hiXt)])f5{y)dy 
{h{y)-E[h{Xt)])fs{y)dy, 



Ae 



■ x{ax + j3) is the squared diffusion coefficient of the Fisher-Snedecor diffusion. 



whereCT2(a;)____ 

Remark 6.1. Solution fh{x) of the equation (|M|) for the identity function h{x) — x has the surprisingly 
simple form: 

h{x) 



1 



a{2- py 

If Y is the random variable with the unknown distribution and Xt ^ TS{a, f3), the Stein equation (|M)) 
takes the form 
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■Y{aY + P)f;^{Y)-9(Y- 



/3 



MY) = h{Y) - E[h{Xt)]. 



a(/3-2) ' '"' ' """^' ' " \' (3-2 

Applying the expectation on the both sides of ([7T|) results in the expression 

29 ,,. _ „. „,.,,. „/,. /3 



(71) 



E 



■ YiaY + /?)/;: (y) -9(y- j^^ MY) 



^ E[h{Y)] - E[h{Xt)]. 



_a(/3-2) ' '"" ' '''"''' ' "V' P-2, 

In particular, according to the general result from Schoutens j64j . random variable Y on (ei, 62) with the 
density function /y(-) and the finite second moment has Fisher-Snedecor distribution if and only if 



E 



29 



a{P~2) 



Y{aY + P)My)-0[Y 



/3 



P-2 



MY) 







(72) 



for all real continuous and piecewise continuously differentiable functions fh on the interval (61,62) for 
which the function g{x) = \2<T'^{x)f'f^{x)\ + \iJ.{x)fh{x)\ is bounded, where /i(x) and 17(2;) are infinitesimal 
parameters of the Fisher-Snedecor diffusion given by (|I1I) . If we define function fh{-) by the expression 

/3-2 



Mx) 



0n{2n - 13) 



f;xx), n = 1, 
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where F'^{x) is the first derivative of the n-th Fisher- Snedecor polynoniial Fn{x), the moment condition 
([7^ takes the form 



/3-2 



9n(2n-/3) 



■E 



29 



a(/3-2) 



Y{aY + li)F':{Y)-e[Y 



/? 



fi-2 



KiY) 



0, n = l,. 



(73) 



According to the Sturm-Liouville differential equation (|A.7p . the final form of the moment condition ([751) 



E[Fr,{Y)]^Q, n=l, 



(74) 



Results of the previous analysis and the following remark will be used for construction of the statistical 
test for Fisher-Snedecor distributional assumptions. 

Remark 6.2. If random variable Y has Fisher-Snedecor distribution with parameters a > 2 and /3 > 4 
then 

E[F.,,{Y)]=Q, n=l,..., ^ 



6.2 Statistical test for known parameters a and j3 

Null hypothesis of the form 

Hq : The probability density function of random variable Xt is fs{x) 

is considered here under assumption that parameters a and /? are known. Furthermore, statistical test 
for testing this distributional assumption about invariant distribution of Fisher-Snedecor diffusion is 
constructed. The basic idea for construction of such a test follows from the GMM approach for testing 
distributional assumptions developed by Bontemps and Maddahi 21i. However, statistical hypothesis 
concerning invariant Fisher-Snedecor distribution is not treated in that paper. 

We consider the sample Xi, . . . , Xn from the ergodic diffusion process {Xt, t > 0} with invariant Fisher- 
Snedecor distribution. Autocorrelation function (fT5|) of Fisher-Snedecor diffusion implies that this sample 
has a decreasing exponential correlation structure. Let us define the function /i: M — >■ M™ which for a; G M 
has a value of the form 



h(x) 



2x{ax + (3)f[ {x) +a[P-{(3- 2)x] h i^) 



2x{ax + (i)f^(x) +a[l3-{(i- 2)x] f^(x) 



(75) 



where /i(-)j ■ ■ ■ i fm{') are real continuous and piecewise continuously differentiable functions. For exam- 
ple, we can define the function h{-) in terms of functions /i(-) — Fi{-), . . . , /m(') — Fm{-), i.e. in terms 
of the Fisher-Snedecor polynomials or even their linear combinations, e.g. 



h{x) = [F^{x), . . . , F^{x)Y, m e {1, . . . , L/3/4J} . 
In this case the Remark 16.21 yields the moment condition 

E[h{Xk)]=Q, Xk^TS{a,l3), 



(76) 



(77) 



which is the basis for construction of the test for Fisher-Snedecor distributional assumptions. The simplest 
case of the moment condition (177)) is based on the function 

h{x)=F,(x), J e {!,..., L/3/4J}, 
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i.e. the moment condition is of the form E[Fj{Xt)] = 0, Xt ^ J-S{a,P). According to the Theorem 15. II 
it follows that 

^^F,(X0 4AA(0,a|), (78) 

where 

E[F,{Xt)]^0, a| =Var(F,(Xo)) + 2^Cov(F,(Xo),F,(Xfc)) = coth(^^^^-^j(/5-2j)j. 

Asymptotic normality (j78p implies that for testing Fisher-Snedecor distributional assumptions we could 
use the following test statistic: 

\-^p^(^^)) -Xl (79) 



'^°*M2(^^^-(/^-2^-) 



Furthermore, we can define the test statistic based on the function (|76p . i.e. in terms of m, e {!,..., [/3/4J}, 
Fisher-Snedecor polynomials. According to the Theorem 15.11 it follows that 

^f]M^t)^AA(O,S(/3,0)), (80) 

where the asymptotic expectation is and the (m x m) covariance matrix 'E{/3,9) — [(jfj(J3,9)]mxm is 
the diagonal matrix whose elements are of the form 

oo oo 

4(/3, 9) = Cov(F,(Xo), F,{Xo}) + ^ Cov{F,{Xo), Fj{Xk)) + ^ Coy{F{Xu), F,{Xo}) = 

k=l fc=l 

= coth {^ ^J^^^ jW - 2j)) S.,, ^,J e {1, . . . , L/3/4J} . 

Then the test statistic for testing Fisher-Snedecor distributional assumptions has the following asymptotic 
distribution: 

xl- (81) 



[-^XH^^)) ^iM-' [-^pi^^)) 



Remark 6.3. Since (3 given by (jBO]) and 9 given by ([56|) are the P-consistent estimators of parameters 
P and 9 and <Jij{(3, 9), i,j £ {1, . . . , [/3/4J}, is continuous function of /3 and 9 it follows that 

[S(^,?)]-i/MS(/3,0)]^/'4l, n^oo. 

Now positive definiteness of the covariance matrix S(/3, 0), asymptotic normality (j80p and Slutsky Lemma 
imply 

n 

Remark 6.4. Proposed goodness of fit test is not consistent against certain fixed alternatives. For 
instance, one can consider an alternative related to the probability density function f{x) with support 
(ei, 62) which is different from fs(x), but such that 

Fnix)f{x)dx = 0, ne{l,...,L/3/4j}. 
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Probability density function f{x) could also be used for construction of such a test. According to the 
nature of our test, it is difficult to expect that it would be consistent against such an alternative. How- 
ever, the goal of this paper is the statistical analysis of Fisher-Snedecor diffusion which is one of six 
Pearson diffusions (see the classification scheme in Introduction). This fact implies that our test is surely 
consistent in this class of stochastic processes. Moreover, since we are dealing with the process with 
known dependence structure it is hard to obtain the distribution of the test statistics that correspond to 
some classical tests. However, these remarks on consistency of our test and its advantages/disadvatages 
in comparison with some classical tests demand further investigation. 



6.3 Testing hypothesis with parameter uncertainty 

The test statistic given in terms of unknown parameters and test statistic given in terms of consistent 
estimators of these parameters can have different asymptotic distributions. Therefore we observe an 
alternative of the statistical test constructed in the Subsection 6.2. in the framework of the parameter 
uncertainty. The approach to this problem is based on the moment condition (j77p , where instead of the 
unknown parameter w = (a,/3) its -^/n-consistent estimator Q = (a,/3) will be used (see Remark [ 
The basic idea is to modify the test statistic ([8T|) by replacing -^ J2^=i ^ i^t) with 



It means that in the parameter uncertainty framework we would like to observe the test statistic 
Asymptotic distributions of the test statistics ([8T|) and ff84|l coincide only if 

^ n ^ n 

-^yh{Xt) and -^yh{Xt,ui) 



(83) 



(84) 



have the same asymptotic distributions. According to the Remark 16.31 it means that we need to check if 

(85) 



-1-1/2 1 

S(^,0) -^Y.^{Xu€:)AN{S),1) 



holds. The classical method for checking the asymptotic normality of (|83p is the Taylor approximation 
around the unknown parameter ui = {a,/3). In particular. 






dh{Xt,Lu) 



y/n{Q -uj) + o,„(l). 



(86) 



where ^ X)"=i — g^V" stands for the sample counterpart of the first moment of the random variable 
If we define 



Mh 



limiV^M^-^=i. 

n— >-oo 77, ^- — ^ 



t=l 



duj^ 



dh{Xt,w) 



the Taylor approximation (1861) could be written as follows: 

'-mxrn Mh 



-^Y.h{X,,€:)^[I„ 



^EtihiX,,^) 
^(cj - uj) 



Om(l), 



(87) 
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where m is the number of the variables of function h and Imxm is the identity matrix. 

Relation ([57)1 implies that the asymptotic distribution of (|55|) depends on the asymptotic distributions 

of 'Te S"=i ^^ i^t,^) — -^ J2t=i ^ i-^t) ^^'^ ^/niu) — uj)~ \/n{a — a, /? — /?) which are given by ((5^ and 

the Theorem I5.2l (iv). respectively. According to Bontemps and Meddahi [21], asymptotic distributions 

of -7=J2"=i^i'^t) and -j=Y^^t^ih{Xt,u)) coincide only if M^ — 0. In that case we say that the test 

statistic (|84|) is robust against the parameter uncertainty. 

To check if this is true, we use the approach based on the generalized information matrix equality (see 

Bontemps and Meddahi [21]). Since the moment condition ((77)) holds for the test function h{-) defined 

in terms of Fisher-Snedecor polynomials, the generalized information matrix equality is 



E 



dh{Xt,io) 



duj'' 



+ E[h^{Xt,^)s{XuUj)] = Mh + E[h^Xt,uj)s{Xt,w)] = 0, 



where s{Xt,ij-!) is the score function of the random variable Xt ^ FS{a,l3). It follows that the condition 
Mh = 0, which ensures the robustness of the test statistic (|84|) against the parameter uncertainty, holds if 
and only if i?[/i'^(Xt, w) s{Xt,uj)] — 0, i.e. if /i'^(Xt,aj) and s(Xf, w) are orthogonal. When Xt ^ TS{a, /S), 
the score function is 



s{x,ui) = 



2{ax+l3] 



2(aa;+/3) 



/3(1 - a;) + {ax + (i) log 



a{x — 1) + [ax + /3) log 



ax+^ 







ax+15 



(ax - 



[ax 



r'(f) r'(^) 




(89) 



From the definition of the test function (|76p in terms of only two Fisher-Snedecor polynomials and 
the expression ([5^ for the score function we conclude that h'^{Xt,uj) and s{Xt,uj) are not orthogonal, 
i.e. Fisher-Snedecor polynomials are not orthogonal to the components of the score function s{x,io). 
It implies that Fisher-Snedecor polynomials are not robust against parameter uncertainty for testing 
hypothesis about Fisher-Snedecor invariant distribution. 
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Appendix A 
Fisher-Snedecor Polynomials 

In this section we refer to some properties of the finite system of polynomials related to the Fisher- 
Snedecor distribution which are called Fisher-Snedecor polynomials. The general theory of orthogonal 
systems of polynomials can be found in the book by Chihara [27] . 

A.l. Pearson equation 

The Pearson equation (see Pearson f6?) for Fisher-Snedecor distribution is the differential equation 

fs'(a;) f3{a- 2) -aiP + 2)x q{x) 



fs(x) 



2x{ax + P) 



s{x)' 



(A.1) 
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Polynomial solutions of the hypergeometric equation 

s{x) y"{x) + [q{x) + s{x)] y (x) + Xy{x) = 0, 

where s{x) and q{x) are polynomials from the Pearson equation (|A.ip . are known as Fisher-Snedecor 
polynomials. 

A. 2. Rodrigues formula 

There exists the unique integer TV G No such that 

oo 

s'^{x)fs{x)dx <oo, ne{0,...,N} and / s"(x) fs(x) dx = oo, ne{N + k,keN}. 

-oo 

In particular, for the Fisher-Snedecor distribution N ~ \j\,(3>2, and for each n e {0, . . . , [/3/4J} the 
function 

Fjx) ^x^-"^ (ax + 13)^+2 J^— (2"xt+"-i(aa; + /3)""^-|\ (A.2) 

ax" I. J 

is the polynomial of at most n-th degree. This formula is known as the Rodrigues formula while the poly- 
nomials defined by (jA.2[) are no n- normalized (with respect to the Fisher-Snedecor density ^) Fisher- 
Snedecor polynomials. The first six polynomials of this type are: 

Foix) = 1, 

Fi{x) ^-a{l3-2)x + al3, 

F2{x) =a^iP- 4)(/3 - 6)x2 - 2a/3(a + 2)(/3 - 4)x + al3^{a + 2), 

%{x) = -a^iP - 6)(^ - 8)(/3 - 10)a;3 + Za'^P{a + 4)(^ - 6)(/3 - %)x^ - ?.al3^{a + 2){a + 4)(/3 - Q)x + 
-Ha/33(a-f 2)(a + 4), 

F^{x) = a^{p - 8)(/3 - 10)(/3 - 12)(/3 - U)x'^ - Aa^^ia + 6)(/3 - 8)(^ - 10)(/3 - 12)a;3 + 
+ 6a2/32(a + 4)(a + 6)(^ - 8)(/3 - 10)a;2 - 4a^3(a + 2)(a + 4)(a + 6)(^ - S)x + 
+ al3^{a + 2){a + 4){a + 6), 

F5(x) = -a5(/3- 10)(/3- 12)(/3- 14)(/3- 16)(/3- 18)a;5 + 5a4/3(a + 8)(^- 10)(/3- 12)(/3- 14)(^- 16)x4 - 
-10a^P^{a+6){a+8){l3-10){f3-l2){P~U)x^+10a^P^{a+4){a+6){a+8){l3-10){P-l2)x^- 
- baP'^ia + 2)(a + 4)(a + 6)(a + 8)(^ - 10)x + a/35(a + 2)(a + 4)(a + 6)(a + 8). 

Fisher-Snedecor polynomials could also be defined in terms of the Gauss hypergeometric function, i.e. 



onon-l t:^ / ^^ -^ /^ . " . " 



i^„(x)=a/3"2"-SFi(-n,n-|;-;--a;l, n € {0, . . . , L/3/4J} . (A.3) 

A. 3. Orthogonality property and normalizing constant 

Fisher-Snedecor polynomials defined by the Rodrigues formula (jA.2p are orthogonal in the general sense 
for n G {0,..., [/3/4J}, /3 > 4. Precisely, there is no set of parameters a and (3 which ensure the 
orthogonality in the restricted sense (see Beale [H]), i.e. 

Fraix)K{x)fsix)dx = 0, m, n £ {0, . . . , [^/4J } , m ^ n. 



To normalize them, we must multiply each of them by the normalizing constant 

(-1)" 
" ^{-irnld^lj 
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where 



dn=Y[ (q'{x) 



fc=0 



n + k + 1 ,,. , 
7; s (x) 



s^\x)fs{x)dx. 



In this particular case d„ and /„ are given by the expressions 



d„ = n [a(2n + 2k - /3)] = (2a)" -^ ^, n G {0, . . . , L/3/4J} . 

fc=o r (^n - f j 



2/32\" i?(n+f,-2n+f 



i3(f4) 



ne{0,...,L/3/4j}, 



while the normalizing constant a„ for Fisher- Snedecor polynomials is given by expression 

a„ = (-1)'' 



r(„-|)i3(|4) 



\ n!(-l)"(2/3)2"rf2n 



2,5 n+f,-2n+f 



(-1)" 



B{jA) 



\ n!(2/3)2"S(f + n, 2 



2n 



.fc=i 



n f+^-2r 



H -1 



(A.4) 



According to (jA.3| and (jA.4p . the Rodrigues formula for the normalized (with respect to the Fisher- 
Snedecor density (j9])) Fisher-Snedecor polynomials is given by the expression 



F^ix) = a„ F„(x) = a„ x^-^ (ax + (3)^+^ — {2" ^^+"-1 {ax + /3)"-t-4 | . (A.5) 

In conclusion, polynomials Fn{x), n e {0, . . . , [/3/4J}, /3 > 4, are normalized and orthogonal, i.e. 



F^{x)Fn{x) fs(x) dx = 5,nn, 



(A.6) 



where 5mn is the standard Kronecker symbol. Relation (jA.6| implies interesting properties of the random 
variables _F„(y), where Y '^ J-S{a, /3). In particular, random variables Fn{Y) are orthonormal, i.e. 

E[FniY)F„,{Y)] - Snm, n, m £ {0, . . . , [l3/4\} . 

Since Fo(x) = 1, substituting n ^ and ?7i = to the previous expression we see that random variables 
Fn{Y) are centered, i.e. 

E[Fn{Y)]=0, ne {!,..., L/3/4J}. 

A.4. Recurrence relation 

Fisher-Snedecor polynomials satisfy the following Favard-Jacobi recurrent relation: 

1 



F„+i{x) = ((6„ - x)Fnix) + a^^iFn-i{x)) , n G {0, . . . , L/3/4J} , 

an 



where 



2(/3-2n) 



a{An - f3){2 + An - f3) 



[n + \){2n + a)(i^{a + /3 - 2n - 2) 
(2 + 4n-/3)(4 + 4n-/3) ' 



/3n(2n -f a - 2) /3(n + l)(2n -h a) 
a(4n - /3 - 2) a(4n - /3 + 2) ' 



fJ-n 



l3{2n + a-2) 
a{An- P-2)' 
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A. 5. Sturm-Liouville equation 

Both non- normalized and normalized, with respect to the density ^, Fisher- Snedecor polynomials satisfy 
the second order self-adjoint equation 

2(ax + B)^^2- d \ x^ ,, A > / N „ 

x^ ' dx y[ax + p)^-^ J 

which could be written as 

2x{ax + P)y"{x) + a (/? - (/3 - 2)x) y'(x) -|- Kv{x) = 0, (A.7) 

where the spectral parameter A„ is given by the expression 

A„ = -n {s'{x) + q{x))' - "^"'^^ ^"(a:) = an (/3 - 2n) , n € {0, . . . , L/3/4J } . 

Fisher-Snedecor polynomials Fi^{x), . . . ,F\ p_\ {x) are the polynomial eigenfunctions of the Sturm-Liouville 
equation (|A.7[) . while real numbers A„ are the corresponding eigenvalues. 

Appendix B 
Gauss Hypergeometric Functions 

Gauss hypergeometric series is the power series 

^ (c)fc fc! c c(c+l) 2 

where z is a complex variable, a, 5 and c are real or complex parameters and {a)k is the Pochhammer 
symbol which denotes the quantity 

(a)o = l, (a)fc = ^-^^-±-^=a(a + l)-...-(a + fc-l), fc G N. 

r(a) 

The series is not defined for c = — ?7i, m ^ Nq, provided that a or 6 is not the negative integer n such that 
n < m. Furthermore, if the series ()90p is defined but a or fe is equal to (— n), n g No, then it terminates 
in a finite number of terms and its sum is then the polynomial of degree n in variable z. Except for this 
case, in which the series is absolutely convergent for \z\ < cx), the radius of absolute convergence of the 
series ([90]) is the unit circle, i.e. \z\ < 1. In this case it is said that the series (|90l) defines the Gauss or 
hypergeometric function 

g{z)= 2Fiia,b;c;z). (91) 

It can be verified (see Slater [68], Section 1.2., page 6) that the function g{z) is the solution of the second 
order differential equation 

z(l - z)g"{z) + {c-{a + b+ l)z) g'{z) ~ abg{z) = 0, (92) 

in the region \z\ < 1. However, the function (PT|) can be analytically continued to the other parts of the 
complex plane, i.e. solutions of the equation ((92)) are also defined outside the unit circle. These solutions 
are provided by following substitutions in the equation ()92|) : 



• substitution z = 1 — y yields solutions valid in the region |1 — zj < 1, 

• substitution z — 1/y yields solutions valid in the region \z\ > 1. 
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Furthermore, solutions in these three regions of the complex plane yield solutions valid in the regions 
|1 — z| > 1, Re(z) > 1/2 and Re(z) < 1/2. The number of solutions in these six regions is 24 and they 
are known as Kummer's system of solutions. This system of solutions is very important since it provides 
the analytic continuation of the hypergeometric function to the whole complex plane with an exception 
of one branch cut which is always the subinterval of the real axis. For example, the analytic continuation 
of the function (|M|) from the region \z\ < 1 to the whole complex plane with an exception of the branch 
cut situated along the interval [l,oo) is provided by the Mellin-Barnes integral (see Slater [55], Section 
1.6.1., Expression 1.6.1.6.): 

zoo 

r(c) f r(« + .)r(6 + .)r(-^^_^^.^^_ ^^3^ 



27rir(a)r(6) J r(c + s) 

— zoo 

Namely, Mellin-Barnes integral defines the analytic function on the complex plane with the branch 
cut along the positive real axis, i.e. along the interval (0,oo). Analytic function (|9ip defined by the 
Gauss hypergeometric series and the analytic function defined by the Mellin-Barnes integral coincide on 
the region {z £ C : \z\ < l,z ^ (0,1)}. Therefore, the Mellin-Barnes integral provides the analytic 
continuation of the function (l?T|) to the region \z\ > 1 cutted along the interval [l,oo). Hence, the 
hypergeometric function 2Fi{a,b]c; z) is an analytic function defined on the complex plane with the 
branch cut situated along the interval [l,oo). Similarly, the hypergeometric function 2Fi{a,b]c; —z) 
is an analytic function defined on the complex plane with the branch cut situated along the interval 
(-00,-1]. 

From the Mellin-Barnes integral follows one of the most important relations for analytic continuation of 
the hypergeometric functions (cf. Whittaker and Watson [73] or Luke [55]): 

g{z) = 2F1 (a, b; c; z) = ^j^j^j^ ~ "j {-zy^ 2^1 f a, 1 - c + a; 1 - 5 + a; i) + (94) 

r(6)r(c-a) V zj 

r(c)r(a-6) , , r, f, , ,1 

''■''■ ' {^z)-''2Fl[b,l-c + b■,l-a + b■,- 



T{a)Tic-b) ^ ' - - \^ ' ■ ' ' '2^ 

where functions (— z)~° 2F1 (a, 1 — c + a; 1 — 6 + a; i) and (— z)^'' 2^1 (fo, 1 — c + 6; 1 — a + 6; i) are solu- 
tions of the equation ([M]) in the region \z\ > 1. A list of useful relations between solutions of equation 
(|92p which providing analytic continuation of the hypergeometric function, can be found in Luke [55) (see 
Section 3.9., page 69). 

Appendix C 
Feller's and Weyl's classification of boundaries of the state space of Fisher-Snedecor 

diffusion 

In this section we present classical Feller's, Weyl's limit-point/limit-circle (LP/LC) and oscillatory/non- 
oscillatory (0/NO) classification of boundaries of the diffusion state space. For general information on 
Feller's classification scheme we refer to Karlin and Taylor [42 and Linetsky [52 , while for Weyl's LP/LC 
and 0/NO classification we refer to Fulton et al. [33 ) and Linetsky [52] • The nature of each boundary 
is determined according to the behavior of the Sturm-Liouville equation (J18p near it. 

Theorem 6.1. Boundaries of the state space of Fisher-Snedecor diffusion with parameters a > 0, a ^ 
{2m, m G N}, and (3 > 2 are classified as follows: 
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(i) Boundary ei — Q is regular for a < 2 and entrance otherwise, while 62 = oo is natural for a > 0, 
a ^ {2m, m g N}. 

(a) For a > 0, a ^ {2m, m G N}, boundary ei — is non-oscillatory, while 62 = 00 is oscilattory/ non- 
oscillatory with unique positive cutoff 

8(/3^2) 
Boundary 62 = 00 is non- oscillatory for A < A and oscillatory for A > A. 

(Hi) Boundary ei = is of limit-circle type for a < 4 and of limit-point type otherwise, while boundary 
62 = 00 is of limit-point type for a > 0, a ^ {2m, m £ N}. 

Proof. The proof consists of three parts, as we observe three boundary classification schemes. 
(i) Feller's boundary classification is based on the behavior of the scale function 

y 
S[x, y] — s(z) dz 

X 

near both boundaries, where 5{x) is the scale density defined by ([T2l) . In case of the Fisher-Snedecor 
diffusion, this non-standard integral is evaluated using Mathematica: 



S[x,y] = 



2-a 



1-— T^ I ^ ^ -, '^ P n '^ Q^ 






2' 2 2' 2' P' 

Since lim 2F1 Tl - f , 1 - f - f ; 2 - f ; -| x") =0 and 2F1 fl - f , 1 - f - f ; 2 - f ; -^ y^ does not 
converge as y — >■ 00, for ot ^ {2m, m G N} the function S\x, y\ has the following properties: 

5(0, y] = lim S\x,y\ < 00, a < 2, 

5(0, y] — lim S[x,y] — 00, a > 2, 
S[x,(X)) — lim S[x,y] — 00, a > 0. 
These results imply that for a ^ {2m, m, e N} and arbitrary e > we have: 

£ e 00 

/ S{0, y]m(y) dy < 00, a < 2, / S{0, y]m(y) dy — 00, a > 2, / S[x, cx))Tn(a;) dx = 00, a > 0, 
£ 

00 € 

/ S[e, y]m(y) dy — cx), a > 0, / S[x, e]m{x) dx < 00, a > 0. 

s 

Therefore, according to the standard Feller's boundary classification scheme, ei = is regular boundary 
for Of < 2 and entrance boundary for a > 2, while 62 = 00 is natural boundary for all positive values of 
a, a ^ {2m, m e N}. 

(ii) Since Ci = is either regular or entrance boundary, it is necessary NO (cf. Linetsky |52|). For 
0/NO classification of the natural boundary 62 = c» the standard procedure is transformation of the 
Sturm-Liouville equation ^TE\\ to the Liouville normal form (cf. Fulton et al. j33j). If we observe the 
Sturm-Liouville equation (jlSp written in the form 

-(xt(aa; + /3)i-t-#/'(x))' = A^^^^a;f-i(ax + /3)-t-|/(x), 
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then 



u{x) = 



^—'■(#)' ^'"'^^-IV^ 



h{u) 






/3cosh ua 



2(/3-2) 



l-a-f! 



— sinh Ua 



2(/3-2) 



Q-l 



and 



QH 



/3^ + (a- l)(a-3)csch^ u. 



[{a + /3f - 1] sech' 



8(/3-2) \^- V ^v ^ \^ ]J2iP-2)^ 

Therefore, the Liouville normal form of the Sturm-Liouville equation (jlSp is 

-5"H + (QH-A)<?(u) = 0. 



2) 



(95) 



In this particular case the natural boundary 62 = oo remains unchanged under the Liouville transforma- 
tion, i.e. the corresponding boundary of the equation (|95|) is 63 = u{e2) — 00. 0/NO classification of the 
boundary 62 = 00 depends on the behavior of the potential function Q{u) near that endpoint. Since 



lim csch ^ 



it follows that 



2(^-2)^ 

lim Q{u) = 



= and lim sech ^ 



2(/?-2) 



= 0, 



:, /5>2. 



8(/3-2)^ 
Since lim Q{u) < 00, according to Fulton et al. (cf. [33j, Theorem 6) 62 = 00 is 0/NO boundary of the 

u— »o o 

equation (|95p with the unique positive cutoff 

flfl2 



A 



:, /3>2, 



8(/?-2)^ 

i.e. it is NO for A < A and O for A > A. According to Linetsky (cf. [S?, Theorem 3), the boundary 
62 = CX3 is non-oscillatory for A = A. This classification of boundaries remains invariant under the 
Liouville transformation (cf. [33], Lemma 2), i.e. 62 = 00 is 0/NO boundary of the Sturm-Liouville 
equation (fT8|) with the unique positive cutoff A. Furthermore, it is NO for A < A and O for A > A. 



(iii) Since ei = is the regular boundary for a < 2, it is necessary of the LC type. For LP/LC 
classification of the boundary ei = for a > 2, a ^ {2m, m G N}, we observe its Liouville trasformation 
el — u{ei) = and liminf u^Q(m) and lim. sup u'^Q{u) (cf. Fulton [33], Theorem 5). Since 



lim u csch u\ 



9 



2(/3-2) 



-^-— — - and lim u sech 
9 u^o ' 



2(13-2) 



= 0, 



it follows that 



and 



Now it follows that 



lim u Q(u) 



liminfw Q{u) =limsupu Q{u) 



(«- 


-l)(a-3) 






4 




,2r) 


r.,N_ ("-!)(«- 


-3) 



liminf M Q(u) < — 
u->o 4 
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for a e (2,4), while 

3 

limsupw Q{u) > — 

u^o 4 

for a > 4, a ^ {2m, m e N}. According to Fulton et al., ei* = is the boundary of LC type for a G (2, 4) 
and of LP type for a > 4, a ^ {2m, m G N}. 

For LP/LC classification of the boundary 62 = oo we observe LP/LC classification of its Liouville trans- 
formation 62 = 00 which depends on the behavior of the function Q{u)/v?' near that endpoint. Since 



lira u ^csch^ u\l —r- — — ■ =0 and hm u ^sech^ u\l —j— — — - = 0, 
M-s-oo \ y 2(p — 2) / ti->oo \ y 2(p — 2) / 



it follows that 



hm^.O. 



According to Fulton et al. (cf. [33], Theorem 7) 62 = 00 is the boundary of LP type for all a > 0, 
a ^ {2?Ti, m g N}. This classification of boundaries remains invariant under the Liouville transformation 
(cf. t33j. Lemma 1), i.e. the boundary ei = is of LC type for a < 4 and of LP type for a > 4, while 
the boundary 62 = c« is of LP type for all a > 0, a ^ {2m, m G N}. D 

Remark 6.5. Since in the Liouville normal form (|95p of the Sturm-Liouville equation (fT5)) the first 
derivative term vanishes, according to Buchholz [24| Wronskian of its linearly independent solutions 

gi{u) = h{u)fi{x(u), -A) = h{u)f4u) and gi{u) = h{u)fi{x(u), -A) = h{u)fi(u) 

is constant. According to the relation ([22|) . it follows that 

gi{u) = h{u) {Bxf3ix{u), -A) + Axf4{x{u), -A)) = B^g^iu) + A^giiu). 

Therefore, the Wronskian of the solutions gi{u) and .94 (m) is given by 

W{g^{u),g,iu)) ^ BxWigi{u),g3{u)) ^ Bxh^{u)WiU{u)j3{u)), 

where B\ is given by the expression (|23p. Expression for h(u) given in the proof of the Theorem 16.11 
implies that 



W(<74(w(x)),5i(u(x)))^2i?;,ai-t/3-^yY^-^^^^ = 2ai-^r^SAAA, a; ^ 00, (C.2) 

where A\ is given by (|2T]) . Since W{g4{u),gi{u)) must be constant, this is its value. On the other hand, 
we know that W{g4{u),gi{u)) = h'^{u)x'{u)W{f4{x{u),~X), fi{x{u),~X)). It implies that 

w^f^ t \ \\ ft f \ \\\ W^(.94(u),5i(u)) 
W{f4,{x{u),~X)Ji{x{u),-X)) = — ^3, , ,, , — 

i.e. 

uifft \\ f( \\\ ^i9i{u{x)),gi{u{x))) 

W{U{x,^Xlh{x,-X)) = ^2(,(,)),,(,(,)) > (C.3) 

h:^(u(x))x'{u{x)) ^ 



where 



s{x) 
and s(x) is the speed density given by (fT3|) . 



Remark 6.6. According to Linetsky (see [52], Lemma 1 and Lemma 2) it follows: 
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(i) for NO boundary ei = there exist the unique (up to the constant factor) non-trivial solution 
(p(x, A) of the Sturm-Liouville equation (1181) which is for all A > square integrable with respect 
to the speed density m{x) in the neighborhood of ei, satisfies the appropriate boundary condition 
at ei and such that (p{x,X) and ip'{x,X) are continuous in x and A in (0,oo)^ and entire in A for 
each fixed x G (0,oo). In our case such a solution is the solution fi{x,~\) given by (|24p . under 
the condition a > 2 which ensures its square integrability in the neighborhood of the boundary 0. 
This solution satisfies the boundary condition 

linr/II:^ = 0. 
x^O s(x) 

(ii) for 0/NO boundary 62 = oo there exist the unique (up to the constant factor) non-trivial solution 
tp{x, A) of the Sturm-Liouville equation p^ which is for all A < A square integrable with respect 
to the speed density m(a;) in the neighborhood of 62, satisfies the appropriate boundary condition 
at 62 and such that ip{x,X) and ip'{x,X) are continuous in x and A in (0, 00) x (0,A) and analytic 
in A < A for each fixed x £ (0,oo). In our case such a solution is the solution f4{x,—X) given by 
(PS)). which satisfies the boundary condition 

x->oo S(X) 
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